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DATA  BASELINE  DEFINITION.  15  JULY  1977 


DATA  ITEM 

SOURCE 

DATE 

1. 

Mass  Properties 

Computer  Listing  PR  *19 

27  January  1977 

2. 

Jet  Control  System 

Rockwell  Block  Diagram  MC476- 

0146.  Rf>v.  F.  Seq.  4 

Received 

29  March  1977 

3. 

Wheel  Control  System 

Rockwell  Block  Diagram  MC476- 

0146,  Rev.  F.  Seq.  4 

Recei  ved 

29  March  1977 

4. 

Earth  Sensor 

Performance  Analysis  of  Static 

Earth  Sensor.  Barnes  Engineering 

Company 

5 August  1975 

5. 

Yaw  Sun  Sensor 

Yaw  Sun  Sensor  Output  Character- 
istics, Rev.  B anti  Rev.  C, 

MCR  CO 56 9 

1 1 February  1 977 

a.  Thermal  Variation- - 

+ 0.6  percent 

(See  Solar  Array  Drive) 

6. 

Thruster  Locations 

Scaled  from  Rockwell  Drawing 

V505-900003,  Rev.  C 

27  January  1975 

7. 

Thruster  Characteristics 

a.  0. 1 lb 

GPS  0.  1 lbF  REA.  Qualification 

Test  Report  CDRL,  Item  A005, 

Contract  F04701 -74-C-0527 

31  March  1976 

b.  5.0  1b 

GPS  5.0  lbF  REA,  Qualification 

Test  Report  CDRL,  Item  A005, 

Contract  F040701 -74-C-0527 

31  March  1976 

i 


LJ  j 


DATA  ITEM 

SOURCE 

DATE 

7. 

c. 

RCS  Tank  Blowdown 

GPS  Propellant  Budget 

Received  from 

Profile,  5.0  lb  and  0.  1 

SAMSO, 

lb  Thruster  Performance 

June  1977 

8. 

Wheel  Dynamics 

Motor  Torque  Speed,  Curve 

14  January  1975 

Plot  No.  5,  Design  Review  Data, 

Item  R-005  Bendix  to  Rockwell 

9. 

Solar  Array  Drive 

3 • 

Pulse  Width  Normal 

Solar  Array  Drive  and  Power 

6 June  1976 

50  + 5 milliseconds 

Transfer  System,  Rockwell 

b. 

Pulse  Width  Extended 

Document  No.  MC476-0148, 

85  + 8.  5 milliseconds 

Rev.  R and  C,  Seq.  05  and  05 

C. 

Slew  Command 

and 

64  pulses 

Global  Positioning  System 

<i. 

Expected  Slew  Normal 

Operations  Handbook,  Volume  1 

11. 85“  + 6° 

Section  2.8:  EPS  Description, 

CDRL  Item  A00G,  Contract  No. 

e. 

Expected  Slew  Extended 

F04701 -74-C-0527 

88.  4"  + 18" 

f. 

Max  Pulse  Rate 

1 to  1. 85 

g. 

Deadband,  1 to  1.7° 

h. 

Driving  Torque 

1 45  oz-in  min 

i . 

Coulomb  Friction 

65  oz-in  max 

j- 

Solar  Array  Drive 

Ball  Brothers  Memo 

5 February  1976 

Torques 

k. 

Minimum  Pulse  Rate 

Phil  Talley  Memo 

22  October  1975 

DATA  ITEM 


SOURCE 


DATE 


10. 

Nutation  Damper 

a.’ 

Radius  of  Curvature 

Blue  Print  Damper  Assembly 

15  January  1976 

158.25"  + 10"  to  18" 

Nutation  GPS  Rockwell  Drawing 

b. 

Length  of  tube  15.75" 

No.  V505-52001 6 

c. 

Mass  of  Ball  .018  slugs 

Rockwell  Viewgraph 

CDR 

d. 

Dissipation  Coefficient 

Honeywell  Memo 

29  September  1975 

c = 0.017  + 30% 

11. 

Pitch  Sun  Sensor 

Ball  Brothers  Test  Data 

Received 

Ball  Brothers  Memo 

15  July  1977 

12. 

Alignment  Specs 

a. 

CES 

Global  Positioning  System  Orbital 

b. 

SSS 

Operations  Handbook,  Volume  1, 

Section  2.4:  Control  System, 

c. 

YSS 

Description  CDRL,  Item  A00G, 

d. 

Reaction  Wheel 

Contract  F04701 -74-C-0527, 

r. 

Thrusters 

P2.  4.  138 

13. 

I’rt 

■cession  Procedure 

a. 

Firing  Pulse  Width 

See  Item  12 

93.  6 msec 

b. 

Expected  Precession 

Angle/Pulse,  0.12° 

r. 

Thrusters  Used 

1 7 and  1 8 

d. 

Firing  Delay,  0 to 

5.24  sec  (increments  of 

1 0.  25  msec) 

e. 

Minimum  Number  of 

Pulses,  4 

DATA  IT  KM 


SOURCE 


DATE 


13.  f.  Procedure:  9 segments 

of  4 minutes  duration; 
reference  sss  pipper 
pulse  each  revolution 

14.  Despin 

a.  Thrusters  Used-- 
11  and  13 

b.  Sequence  Modified  from  Item  12 

95  to  85  rpm 

Calibrate 
85  to  75  rpm 

Recalibrate 
75  to  20  rpm 
20  to  10  rpm 
10  to  6 rpm 
6 to  1 rpm 

c.  Firing  Time  in 
Increments  of  5.  27  sec 

15.  Earth  Acquisition  Procedure 
at  1 rpm 

a.  Enable  roll  axis  control  See  Item  12 

b.  Wait  5.  0 minutes 

c.  Enable  pitch  axis  control 

16.  -y  Panel  Deployment  See  Item  12 

Procedure  with  roll  and  pitch 

jet  control,  deploy  -y  panel, 
slew  90  deg 


SAMSO  Review 
15  July  1977 
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DATA  ITEM 


SOURCE 


DATE 


17.  Sun  Acquisition  See  Item  12 

At  0.  75  rpm  and  -y  panel 

deployed  and  sun  line  parallel 
to  panel  normal,  enable  yaw 
jet  control. 

18.  +y  Panel  Deployment  with  See  Item  12 
three-axis  jet  control, 

deploy  +y  panel,  slew 
90  deg 

19.  Momentum  Dump  See  Item  12 

Procedure 

a.  Automatic  dump  of  pitch 
and  roll  axis  when 
momentum  is  greater 
than  0.  25  ft-lb-sec 

b.  Pitch  axis  firing  times 
0.1,  0. 3,  0. 5,  or  0.  7 
sec  (0.3  nominal) 

c. 


Roll  axis  firing  times 
0.8,  2.  4,  4.  0,  or  5.  6 
sec  (2.4  nominal) 


APPENDIX  B 


ONE- BODY  MODEL  EQUATIONS 


This  appendix  presents  the  derivation  of  equations  of  motion  used  in  the 
one  body  simulation.  For  Definition  of  Symbols,  see  Figure  B-l. 


Vector  torque  equations  are  given  by 
/ rxdF  = / rx  dm 


- d2r 
rx  — » 

dt 


d - - , dr. 

= — - / rx  (-jr)  dm 
at  dt 


(B-l) 


t = 

1 dt 


where  H is  defined  as 


H = / rx  (-^)  dm 


= [J  + m M(r  ,r  ) + Z J.  + E m M(r  ,r. )]  u + E J.w, 
L o o oo  ^ 1 ^i  11  o ^ii 


- I a + El . uj . 
o o i 1 1 


(B-2) 


where , in  body  0 coordinates 


AT  T 

M(a  ,b)  = a bl  - a b 


-I 


-I 


xy 


xz 


xy 


xy 


-I 


xz 


xy 


J = T J (T  is  transformation  from  wheel  i to  body  o) 

i oi  i io  oi 
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o) . — T ,oj  . 
i oi  1 


TX  - f (y2  + z2)dra  = Xyy  + 1ZZ  * PtC- 


I = / xy  dm,  etc. 
Ay 


The  rotational  equations  of  motion  are  then  given  by 


T = 


dH 

dt 


= -rr  ri  oi  + z i « ] 

dt  L o o . i iJ 


5 t 


[lo“o  * \ W + “o*  [lo"o  + l h“l] 


T = Vo  + l Vi  + “o3* 


Wheel  Torques  may  be  described  by 


/ x.xdF  = / x.  x — r x dm 
dt 


d dx^ 

= dF  / xix  dT  dm 


dH. 


T = — — - 
x dt 


Wheel  momentum  is  given  by 


H . = / x.x  — x.  dm 
x i dt  l 


= / xlX  x.  + x xA] 


dm 


10 


J 


f x^x[w^  x x^  + wq  x x^]  dm 


= J.co  + J tu 
i i i o 


(B-4) 


Differentiating 

T,  = + J,<I>  + u>  x J.w.  + u x J.oj 

i ii  io  o ii  o 10 


(B-5) 


Torques  about  the  wheel  axis  are  given  by 


T.  = J.w.  + J.i 
1 l i i o 


(B-6) 


Vi 


(B-7) 


where  T^  is  the  sum  of  motor,  bearing  and  windage  torques  on  wheel  i. 
Equations  B-7  and  B-3  are  combined  to  compute  the  body  rate  derivatives. 


Id  + t T + <ju  xH 
o o ^ i o 


(B-8) 


Kinematics 


The  inertial  coordinate  system  and  the  Euler  sequence  are  consistent  with  the 
Rockwell  convention  presented  in  the  AVCS  splinter  meeting  notes  at  the 
preliminary  design  review.  Specifically,  Rockwell  specified  a yaw  pitch 
roll  Euler  sequence  with  the  equations: 

0 = u>  cos  y + q cos  (J>  - r sin  <t> 

ij,  = (w  sin  ip  sin  6 + q sin  <p  + r cos  <p)  (1/cos  6) 

• • 

t = p + ij)  sin  6 + in  sin  i p cos  6 
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whe  re 


to  = orbit  rate 
P >r  = body  rates 


1 


<t> , Q , = Euler  angles 
4> , 6 , ^ = Euler  rates 

These  equations  are  consistent  with : 

1.  A local  vertical  system  with  z axis  pointing  toward  the  earth  and  x 
axis  in  the  orbit  plane  and  positive  on  the  direction  of  motion. 

2.  Euler  angles  specifying  body  orientation  relative  to  the  local  vertical 
frame . 

The  inertial  coordinate  system  is  the  Geocentric  Equatorial  Frame  (e.g.,  z 
axis  positive  North  and  x axis  positive  in  the  direction  of  the  Vernal 
Equinox).  The  local  vertical  system  is  consistent  with  Rockwell's  (e.g., 
z axis  positive  toward  the  earth  and  x axis  in  the  orbit  plane). 

Two  sets  of  Euler  angles,  both  having  a yaw  pitch  roll  sequence,  are  used. 

The  first  set  PSI  , THETA1,  PHI1  define  the  inertial  to  body  transformation  E. 
The  second  set,  PSILV,  THETA,  PHI  specify  the  orientation  of  the  body 
relative  to  the  local  vertical  frame  and  should  be  consistent  with  Rockwell's 
Euler  angles.  The  inertial  Euler  angles  are  integrated  in  the  simulation 
via  quaternions. 

Linear  Motion 


The  simulation  uses  two  subroutines ( POSITN  and  RADVEC)  to  compute  the 
ephemerides  for  the  satellite  and  the  sun.  The  ephemerides  for  the  satellite 
and  the  sun  are  calculated  using  the  six  Kepler  elements:  semi  major  axis  (a), 


* 
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eccentricity  (e)  , inclination  (i)  , right  ascension  of  the  ascending  node  (U)  , 
argument  of  perigee  (w) , and  mean  anomaly  (M) . The  update  is  a simple  rate 
term  on  the  latter  three.  Table  B-l  presents  the  six  elements,  their  units 
and  values  for  the  vernal  equinox  condition.  The  rate  terms  are  derived 
in  the  program  from  the  fixed  parameters  and  are  given  by 


k"\/  • /a  ^ 2 

u = '■  # "9-  (2.5  sinZ  i-2) 

[a(l-e2)]2 

where 

p = Gm  = 62630.3949 

2 

k = Rg  J 2 (first  term  of  perturbing  force  = -19255.96124 
due  to  Earth's  oblateness) 


The  rate  terms  for  the  sun  are  fixed  and  are  included  in  Table  B-l.  The 
simulation  updates  M,  fi,  u by 

M(t)  = M(t  ) + M • (t-t  ) - - < M < T. 

o o — 

ft(t)  = U(.t  ) + h ’ (t-t  ) — TT  < fi  < - 

0 0 — 

0)(t)  = U)(t  ) + U)  (t-t  ) — IT  < U)  < IT 

O O — 

A seventh  term,  the  Greenwich  Hour  Angle  (GHA)  , which  relates  the  earth- 
fixed  coordinate  system  to  the  inertial  system,  given  the  six  elements  is 
updated  via: 


GHA(t)  = GHA( t ) + w ‘ (t-t  ) 
o e o 


a!  = earth  rotation  rate 
e 
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TABLE  B-l.  GPS  AND  SOLAR  EPHEMERIDES 


STANDARI 

) VALUE 

SYMBOL 

NAME 

UNITS 

GPS 

SUN 

a 

Send  Major  Axis 

Miles 

14342. 

1* 

e 

Eccentri ci ty 

— 

0. 

.0167 

i 

Inclination 

Dc  grees 

63. 

23.443 

M 

Mean  Anomaly 

Degrees 

0. 

88.068 

ft 

Right  ascension 

of  ascending 

node 

Degrees 

0. 

0. 

0) 

Argument  of 

perigee 

Degrees 

0. 

270. 

-5 

M 

Degrees /sec 

— 

1.14074x10 

ft 

Degrees /sec 

— 

. 544  773xlO~9 

CO 

Degrees /sec 

— 

0.0 

GHA 

Greenwich  Hour 

Angle 

i 

De grees 

0.0 

Astronomical  unit 


APPENDIX  C 

THREE-BODY  DYNAMIC  MODEL 


SOLAR  PANEL  PIVOT  POINTS  Pj 
BODY  CM's.  CIVIj,  CMj,  CMk;  CM  TOTAL  CM 
MASSES:  nij,  mj,  mk,  mn;  m = TOTAL  MASS 
INERTIA  TENSORS  (ABOUT  CMs)  J„,  Jj,  Jj.  •>!<.  >*„ 


Figure  C-l.  Notation 


+ v I*  (r  + r + a + x ) x + IL  x (r  + r.  + p.  + x.)  + o>.  x (p.  + x ) dm 

^ j ° j J j 1 fit  o o ] 3 3 3 3 / 


6(r  +p  + x ) 
onn 

£ ; (r  + r +p  +x)x 

J o n n n fit 

n n L 


+ a x (r  +r  + p +x  ) dm 
o o n n n 


Noting  that 


m ? + £ m. (r  + r.)  + £m.(r  + r.  + p ) + £m  (r  + r + p ) = 0 

oo.ioijjo3  3nno  n n 


Equation  (C-2)  becomes 


11  = I in  + £ I.u>.  + £ I. to.  + £ H 
o o .11  . J n 

l j J n 


where 


I = J + £ J + £ J.  + £ J - mM(r  , r ) + £m.M(r.,  r.) 
oo.i.]  n 00.111 

i j J n i 

+ £m.M(r.  + p.,r.  + Pj)  + + P„.  ^ + P„> 


I.  = J. 

l l 


I.  = J.  + m.M(r  + r.  + p.,  p.) 

3 J 3 ° J J 3 


H = m (r  + r + p ) x (fip  /fi()  + f (r  + r + p + x ) x (fix  /fi  ) 
n n o n Mn  n t J o r n n n t 


AT  . T 
M(a,b)  = b al  - ba 


COMPUTATION  OF  II 


The  vectors  ? , ? , r , are  given.  The  vector  r can  be  computed  from  Equation 
i 3 n I ° 

(C-3). 


Wheel  Momenta 


I.u).  - T .j!  u). 
11  oi  1 i 


Panel  Momenta 


I.iD.  = T . J.  (i) . + m.M(r  + r + p , p ) m 
3 3 03  3 3 3 o i 3 3 3 
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where 


0 


t K 


U).  = T . / 0 . >=  T . ctf 
3 oA  p]  ( 03  3 


N utation  Dampers 


Assume:  J = 0 (point  mass) 

" 6P, 


H = m (?  + r + p ) x 
n n o n n fit 


ROTATIONAL  DIFFERENTIAL  EQUATION  IN  SATELLITE  COORDINATES 


— dll  SH  - 5 

T = — = — — +11)  X H 
dt  6t  o 


T = I to  + i cu  + El.i.  + £ Li).  + £1.0).  + £ H 
00  00  j11  j 3 3 j 3 J n" 


EVALUATION  OF  THE  DERIVATIVES 


Inertia  Tensors 


• m • T 

I =7.3.-  m[M(r  , r ) + M (r  , r )]  + Em  [M(r  + a,  p.)  + M (r,  + P-.  ft>l 

o.j  OO  OO  • J J J J J J J 

1 J 

+ E mn[M(rn  + P„.  P„)  + MT<rn  + pn>  pn)] 


i =3.  +m.M(r  + r.  + p.,  p.)  + m M(r  + p , p.) 
J J J o 3 3 3 3 0 3 J 


where 


p.  = 0).  x p. 
3 3 3 


ro  m 


Im.p.  + Em  p 
7 33  n n 

3 n 


WHEEL  TORQUES 


J x.  x (IF  - f x.  x — j x 
. 1 v 1 jt2  1 

1 1 dt 


T.  = J.iu  + J.ii . + u>  x [J.(o)  + u) .)] 
1 10  11  O to  l 
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1 


SM. 


. j . .a......  . . 


The  magnitude  of  the  wheel  torques  about  Ihe  axis  is  given  by 


T.  = h.  • T. 
ia  i i 


where  h.  is  a unit  vector  along  the  spin  axis.  For  example. 


h.  = T . 
1 oi 


PANEL  TORQUES 


I (,1j  + xj>  x <IF  = J <Pj  + xj)  x ^ <Pj  + x.  + rp.)  dm 


d2r 

T.  = mp-.x-il  +±Lj 

' J dt2  (,t 


where 

H = f.T  + m M(p  , p )]  (d  + i ) 
j J In  .1  1 1 o 

dr  . 

T . = m . p.  x — J2A 


j ’ "VP  + tJj 4 nyMlcr  v 4 « «,m  (.j  * »„> 

+ [J.  + m M(p  , p .)]  (<».  + i)  ) + u>  x H 
J J J J jo  ° J 


r = h.  • t. 

J j 


where 

h.  = unit  vector  along  panel  hinge 


= T . 
oj 


and 


V 

d r . 

P-  = p x a 

+ p . X 

[A 

dt2  -*  cm 

j 

L «2 

+ M(p.,  r 

+ r.)  i 

+ a)  : 

J o 

J o 

o 

8 r 

— ^ = ‘ m 

P.  + £ m o 1 

5t“  j 3 

J m 

n n 

6 r 


+ tu  x 


+ [M(p.,r  ) + M(p.,r 
jo  jo 
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APPENDIX  D 

NUTATION  DAMPER  MODEL 
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APPENDIX  D 

NUTATION  DAMPER  MODEL 

This  appendix  summarizes  the  nutation  damper  model  implemented  for  the 
digital  simulation  evaluation  of  GPS  stability  and  control  analysis. 

The  nutation  damper  is  used  for  passive  GPS  spin  axis  stabilization  during 
the  initial  mission  phases  by  dissipating  the  spacecraft  precession  energy. 
The  nutation  damper  adopted  for  GPS  application  is  of  the  type  of  ball  in 
gas  filled  tube.  The  tube  is  made  of  aluminum  and  the  ball  of  tungsten 
carbide.  The  gas  is  a 9 : 1 ratio  mixture  of  nitrogen  and  helium  at  1 atm. 
pressure.  Physical  parameters  of  this  danper  is  shown  in  Figure  D-l. 

A similar  damper  of  this  type  was  used  in  the  Telstar  satellite,  Reference 
. D-l.  Discussions  on  the  design  and  analysis  of  ball-in-tube  type  nutation 

• dampers  can  be  found  in  Reference  D-2. 

Derivation  of  the  torque  equations  presented  here  follows  the  same  procedure 
of  Reference  D-2  which  contains  numerous  typographical  errors.  The  results, 
in  the  torque  form,  are  compatible  for  integration  into  the  three  body  model 
described  in  Appendix  C. 

I 

DAMPER  MODEL  DERIVATION 

To  facilitate  the  model  derivation,  we  shall  define  a vehicle  fixed  B-frame 
whose  origin  is  at  the  vehicle's  center  of  mass  (excluding  damper  ball)  with 
the  (X  - Y - Z ) axes  aligned  with  its  principal  axes.  Also  in  the  fol- 

— D — D — D 

lowing  consideration  of  vehicle  rotational  motions,  it  is  assumed  that  the 
translational  motion  of  vehicle's  center  of  mass  resulted  from  reactional 


IN  fV\i  FILLED  TUBE 
D TG  PROVIDE  CcNTRI  RIGA  SPRING 


Figure  D-l.  Nutation  Damper  Characteristics 


forces  of  the  danger  ball  can  be  neglected.  By  reference  to  the  geometrical 
description  of  damper  location  in  B- frame,  the  position  of  center  of  ball, 

R,  can  be  written  as.  Figure  D-2, 


<iV  ■ 


-a  + Y cos  a 
b + Y sin  a 


(D-l) 


whe  re 

Y = 159  in.  = 13.25  ft. 

(a,b)  = coordinates  of  tube  center  of  curvature 
Equation  (1)  inplies  that  Y - 7.  plane  contains  the  damper  tube, 

D D 

Let  the  vehicle  moment  of  inertia  be  such  that 
PxX  = *YY  = A 


(D-2) 


The  rotational  energy  of  the  vehicle  becomes, 
Tv  = % A(WX2  + Wy2)  + %CWZ2 


CD-3) 


The  velocity  of  damper  ball  can  be  written  through  differentiation  of  (1)  as 

b)  = TBI(5  b)  = TBI^-  b^ 


with  (RIb)  = TfiIT  (RBb)  , we  have 

= tbiT  <?b>  + tbiT 


= -trtT  s [vW8,)]  (RBk)  + tbit  (lBb) 


ZB  = SPIN  AXIS 

cMv  = VEHICLE  CENTER 
OF  MASS 

DAMPER  CONTAINED  IN 
PLANE  (Yr  - ZR) 


Figure  D-2.  Nutation  Damper  in  Vehicle 
Body  Fixed  Frame 


Cf\)--sr(5\)]  (RB. ) + (t\> 


In  component  forms,  (4)  becomes, 


o -wz  WY 

wz  0 “wx 

-wY  wx  0 


-a  + y cos  a 
b + y sin  a 


+ — y sin  ct  a 

Y cos  a a 


Hence  the  translational  kinetic  energy  of  the  damper  ball  is 


Tb  - % M <VX2  + Vy2  + Vz2) 


= H M { £ — ( — a + y cos  a)  + (b  + 
+ [-(b  + y sin  a)  Wy  ~ Y s:*-n 


a)  Wz  + (b  + y sin  a)  Wy]' 


•l2 

a aj 


(D-4) 


(D-6) 


y 

+ [(-a  + y cos  a)  wx  + Y cos  a & } 

where  M = mass  of  the  ball 

= 7TYb^  p = 26  3 gm  = 0.018  slug 

3 

p = density  of  tungsten  carbide  (WC)  = 15.7  gm/cm 


Assuming  the  ball  is  rolling  without  slipping  in  the  tube,  the  magnitude  of 

the  angular  velocity  of  the  ball  is  related  to  the  translational  velocity 

of  the  center  of  the  ball  as 

„ 1 r^B  'J  /5b 

Wb  ~^b}  ^b^ 

1 _ 2 , 2 .2,2  2 .2-h 

- rySinaa^Y  cos  a a J 

^b 


= -*-a  where  Y,  = radius  of  the  ball  = 0.625  in. 
>b  b 


(D-7) 


28 


Hence,  the  rotational  energy  of  the  ball  is 
T , 4 <4  M Yb2)  d±>2 

= J M Y2  a2 


(D-8) 


The  drag  on  the  ball  is  viscous.  Reference  D-3,  since  the  tube  is  gas  filled, 
the  clearance  between  ball  and  tube  is  small  and  the  ball  moves  much  slower 
than  the  speed  of  sound.  Hence,  the  viscous  drag  on  the  ball  can  be 
represented  by  Rayleigh's  dissipation  function  as 


f ■ >*D  <iVT  <sV 

, „ 2 .2 

= dCp  y a 

where  is  a constant  to  be  described  later. 


(D-9) 


The  total  energy  of  the  system  then  becomes 
T = T + T,  + T 

v b y 

? 2 2 2 

= W {4A  + (y  + b +a  + 2b  y sin  a -2a  y cos  a)} 

X 

2 2 2 2 

+Wy  {4A  + (b  + 2b  Y sin  a + y sin  a)} 

2 2 2 2 

+WZ  {4C  + 4M  (a  - 2a  y cos  a + y cos  a)} 

2 

+WYWZ  (_ab  " a Y sin  a + b Y cos  a + y sin  a cos 


+W  a (M  (y  -ay  cos  a + b y sin  a) } 

A 


+ -2  M 

+a  {10  My  } 


(D-10) 


With  the  total  vehicle  energy  given  in  (P-10)  and  the  damper  dissipation  func- 
tion given  in  (D-d),  the  equation  of  motion  of  the  vehicle  and  the  damper 
system  can  be  obtained  as 
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The  rotational  equation  of  motion  of  the  vehicle  can  be  written  in  the  form 
of  Euler's  rigid  body  equation  as 


AWX  + (C  - A)  WzWy  = Tx 
AWy  - (C  - A)  WXWZ  = Ty 

AWZ  = Tz  (D-12) 

where  the  torques  Ty,  Ty , and  Tz  are, 

= -M  [Wx  (y^  + b^  + + 2b  y sin  a - 2a  y cos  a) 

+a  (2W  + a)  (by  cos  a + a y sin  .) 

A 

2 

+a  (y  -ay  cos  a + b y sin  a) 

2 2 2 2 2 2 
+WyW  (a  - b - 2a  y cos  a - 2b  y sin  a - y sin  a + y cos  a) 

2 2 2 
+ (Wy  - W ) (ab  + a y sin  a - b y cos  a - y sin  a cos  a )]  (D-13) 


2 2 2 2 
+ W W (y  + b + 2b  y sin  a - y cos  a) 
X 

i 2 

+ 2W^  a (b  y cos  a + Y sin  a cos  a) 

2 2 

+ 2WZ  a (b  y sin  a + Y sin  a)] 

2 2 2 

= -M  [Wz  (a^  - 2a  y cos  a + y cos  a) 

+ (Wv  + W W ) (ab  + a y sin  a - b y cos  a 
l X X 

2 2 2 
+ W Wv  (-a  + 2a  y cos  a - y cos  a) 

x i 

2 

+ 2WZ  a (a  y sin  a - y sin  a cos  a) 

2 2 

- 2W^  a (b  y sic  a + Y sin  a) ] 


(D-14) 


Y sin  a cos  a) 


(D-15) 


With  ball's  equation  of  motion  in  tube  given  as, 

7 . 

— Ya  + W (y  ~ a cos  a + b sm  a) 

D X 

2 2 

- W (b  cos  a + a sin  a)  - W (b  sin  a + y sin  a cos  a) 

X ^ 

2 

- Wz  (a  sin  a - y sin  a cos  a) 

2 2 

- WyWz  (a  cos  a+b  sin  a-  y cos  a + y sin  a) 


+ W a (a  sin  a + b cos  a - a sin  a - b cos  a) 
X 


+ (ir  > « = 0 


(D-16) 


The  drag  force  on  the  ball  consists  of  two  components:  the  pressure  drag 
and  the  viscous  drag.  For  the  case  where  the  gap  between  tube  and  ball  is 
small,  as  in  the  GPS  nutation  damper,  and  the  speed  of  the  ball  in  the  tube  is 
small  conpared  with  the  speed  of  sound,  it  has  been  determined  in  Ref.  D-3 
that  the  drag  force  is  primarily  due  to  pressure  drag.  The  drag  force  D, 
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T 


can  be  written  as 

/ 3 


D = 


135 

4x64 


where 

u = viscosity  of  gas  in  tube 

g = (ya  - Yb>Yb  = gap  parameter 

v.  = inside  radius  of  tube 
' l 

S = speed  of  ball  relative  to  gas 
0^  = dissipation  coefficient 

135  ^ p Yb 

4x64  5/2 

g 


For  GPS  nutation  damper, 

p = 178.1  x 10-6  x 2.089  x 10~ 3 lbf  sec/ft2 
g = (0.633  - 0.625)  / 0.625 
y . = 0.633  in. 


Hence,  value  of  dissipation  coefficient, 

CD  = .0171  lb  f /(  ft  /sec) 

To  handle  the  discontinuity  at  either  end  of  the  damper  tube,  let 
d = -ed  at  a = + a max 

where 

a max  s (8/159)  rad  =2.9  deg. 
e =0  for  inelastic  collision 
= 1 for  elastic  collision. 
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The  GPS  requirement  on  half  cone  angle  of  nutation  is  to  be  less  than 
2 deg.  The  discontinuity  situation  is  not  expected  to  be  encountered  often 
if  the  2 deg.  requirement  is  met.  The  value  of  e = 1 would  give  a 
conservative  result  in  that  less  precession  energy  is  dissipated  at  bottoming. 
Hence,  use  of  c = 1 is  recommended  for  this  reason. 

NUTATION  DAMPING  TIME  CONSTANT  VERIFICATION 

Nominal  Solution  Undamped  Rigid  Body  Motion 


Consider  a rigid  body  with  symmetric  moment  of  inertia  about  the  traverse 
axes,  i.e.  , I = 1^  in  a torque  free  situation.  The  Euler's  equation 
is  reduced  to: 


The  angular  momentum  vector  is  written  as: 


- ' ZXXWX  -B  + IYY  WY  -B  + IZZ  WZ  h, 


(D-20) 


Since  there  is  no  external  torque,  the  angular  momentum  is  constant.  The 
nominal  motion  of  the  rigid  body  can  be  visualized  as  a constant  precession 
about  the  angular  momentum  vector  with  a precession  rate  Si  and  a nutation 
angle  0 which  is  defined  as: 


tan  0 = — = 
HS 


2 2 2 
+ I,™  W z 
YY  Y 


Tzzwz 


iyyb 


(D-21) 


Therefore ; 


77 

~ x tan  9 I W 
-**yy  “ 


For  small  value  of  0,  (6)  becomes 
I lr 


(D-22) 


B 


ZZ 

CYY 


x 6 


W„ 


(D-23) 


A relationship  between  the  rigid  body's  angular  momentum  and  rotational  energy 
useful  for  later  development  of  the  energy  sink  approximation  is: 


2 2 

9 ? o Hf-  H 

T = hi  (WZ  + W)  + m Wz  = — H 2_ 

2 XX  V X Y ; ‘ ZZ  Z 2I„„  21 


(D-24) 


XX 


ZZ 


2 2 2 
H = H + H 
t s 


Hence  H. 


From  ( 5) 


I I 

/ 2 1 

XX  ZZ 

2T  - H— 

I - I 

z I 

ZZ  XX 

ZZ 

(D-25) 


. 2 Ht 

sm  ■ = — y 
H 


IxxIzz  ' 


ZZ 


xx 


2T 

„2 


ZZ 


(D-26) 
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_ M 


Nutation  Damping  Time  Constant 


The  GPS  nutation  damper  consists  of  a tube  aligned  with  respect  to  the 
vehicle  spin  axis.  A ball  inside  the  tube  is  dissipating  the  nutation 
energy  via  viscous  drag  force  retarding  the  ball's  notion.  An  approximation 
procedure  known  as  the  energy  sink  concept  assumes  that  the  dissipation  of 
nutation  energy  does  not  change  the  system's  angular  momentum  (which  is 
dominantly  spin  angular  momentum) . Hence  equation  (0-2 A)  car  be  differentiate 
with  terms  H neglected  in  this  approximation,  i.e., 

21  T 

2 sin  0 cos  0 6 = — ~ (D-27) 

(n-l)H 


For  small  value  of  0, 


e 0 


2I„„T 

ZZ 

(n-l)F' 


A further  approximation  involves 


2 2 
H " HS 


2 2 
I W 
ZZ  Z 


(D-28) 


(D-29) 


Hence 


<"-1>Izz“z  0 


(n-30) 


The  energy  dissipation  T can  be  expressed  as  the  product  of  the  drag  force 
the  ball  experienced  and  the  velocity  of  the  ball,  i.e., 

T = -fv  (0-31) 
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with  drag  force, 

f = cDv 

and  the  ball  velocity. 


V = ya 


Hence  the  energy  dissipation  rate, 

• 2 • 2 

T = -CD  ' 01 


(D-32) 


where  a is  the  position  of  ball  in  tube  defined  as  in  Figure  D-2, 

Y = is  radius  of  damper  tube, 

CD  = is  dissipation  coefficient 

From  analysis  of  nutation  danper,  the  motion  of  the  ball  in  tube  was 
obtained  as 

7 2 

yot  + Wx(y  - a cos  a + b sin  a)  - Wy  (b  cos  a + y sin  a cos  a) 

2 2 

-W  (b  cos  a + a sin  a)  - (a  sin  a - y sin  a cos  a) 

2 2 

-l^W^Ca  cos  a + b sin  a - y cos  a + y sin  a) 

CD7 

+ -jjT  “ = 0 ( D-  3 3 ) 


With  the  geometry  of  mounting  such  that  b = 0 and  introducing  the  small 

2 2 

angle  approximations  sin  a ~ a,  cos  a = 1 and  neglecting  terms  W and  Wv 
as  they  are  small  for  small  nutation  angle  O,  equation  (17)  becomes 


(t^)«  + ("TrOa  + W„2(y  - a)a  = WVW  (a  - y)  - W (y  - a) 


(D-34) 


Using  the  expression  of  Wy  as  given  in  equations  (2)  and  (7) 


on  Wz  cos  ut 


(D-35) 


Similarly , 


= -6n  fi  cos  fit 


Equation  (18)  can  be  reduced  to  the  following  form: 


. 5Sl  • f5<Y  - a)V 

"MaT  5+  1: F J“ 

5(y  - a)  (2n  - n2)W  2 

- t T~< 


cos  fit 


(D-36) 


The  solution  of  the  standard  second  order  linear  differential  equation  of 


the  form 


+ 2b  a + a ct  = c9  cos  fit 


(D-37) 


has  been  obtained  in  standard  textbooks  as: 


, — 5 ~ cos  (fit  +<?>)  = 6a  cos  (fit  + $) 

(a  - fi  ) + 4b  fi 


(D-38) 


Compare  (20)  and  (21) , 


2b  = 5 Cd/(7M) 

a2  = 5 (y  - a)Wz2/(7y) 

c = -5 (y  - a)  (2n  - n2)W  2/(7y) 


CD-39) 


From  (22) , 


a = - 8 a fi  sin  (fit  + $) 
o 


Substitute  (24)  into  (16) 

. ? 9 9 9 2 

T = -0  C v u fi  sin  (fit  +i) 
D o 


(D-40) 


(D-Al) 


The  average  value  over  a full  period,  i.e.. 


2 e2 


fit  = 2*  , is 

P 

T = -h  c Y2  a 2 Q2  02 

D o 


1 9 

since  — / P sin2  (fit  + 0) dt  = h 
P 

Substitution  of  (26)  into  (14)  gives 

® — 1 Cp  2 a 4 

l_(n-l)lz2wz2ej  L 


2 2 2 

S'-  °o  ° 

i(n-mzz»z2 


The  time  constant  of  nutation  damping  can  now  be  written 

T 2fo-1>Izz"z2 
c S-'2-2'-'2 

D o 

where 


(D-42) 


(D-43) 


(D-44) 


"°  ' V (,2  - 22)2  + 4b2  !,2 


2b  = 5 Cd/(7M) 

a2  = 5(y  - a)Wz2/(7y) 

c = -(2n  - n2)a2 


fi  = (n  - 1)W„ 


The  predicted  nutation  damping  time  constant  for  the  two  cases  of  Wr 
100  RPM  and  W = 10  RPM  are  computed  as  follows: 

Case  1 : 

r = 13.25  ft. 

a = 10.25  ft. 

m = 0.018  slug 

CD  = 0.0171  lb  f/(ft/sec) 

1XX  = XYY  = 77  slu^ft2 

TZZ  = slu8”ft2 

Wz  = 100  RPM  = 10.47  rad/sec. 

n = 1.22 
SI  =2.31 

2b  = 0.6  789 

a2  = 17.728 

c2  = -16.84 

a = -1.  348 
o 

T = 155.755  sec. 

c 

Case  2 : 

Wz  = 10  RPM  = 1.047  rad/sec. 
n = 1.22 

£2  = 0 . 2 31 

2b  =0.6786 

2 
a 


0.17728 


c 


= -0.1684 


a = -0.84275 
o 

T = 398.497  sec. 
c 

The  nutation  half  cone  angle  history  corresponding  to  Case  1 and  Case  2 
considered  above  are  plotted  in  Figures  D-3  and  d-4  respectively.  These 
results  are  obtained  using  die  nutation  damping  simulation  program  where  the 
detailed  dynamics  of  the  damper  ball  and  rigid  body  dynamics  are 
consi dered . 


4C 


i 

> 


100  RPM 


*xx  = 76-3 

Jyy  = 78-3 

1 zz  = 94-1 


"Time  constant"  5 150  sec. 

(155.75  sec.  was  predicted) 


Figure  D-3.  Simulation  Data 
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10  RPM 


xx 

'yy 


ZZ 


76.3 

78.3 

94.1 


M 


Time  Constant"  * 365  sec. 

( 398.5  sec.  was  predicted) 
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APPENDIX  E 


MODEL  FOR  GPS  DIGITAL  SPIN  SUN  SENSOR  (SSS) 


This  appendix  documents  the  GPS  spin  sun  sensor  (SSS)  model.  It  has  been  programmed 
on  a digital  computer,  and  the  simulation  performs  satisfactorily.  Figure  E-l  shows 
a flowchart  of  the  math  model  which  is  based  on  the  data  of  the  references  and  the 
assumptions  listed  in  this  appendix.  The  main  simulation  program  will  provide  the 
model  with  the  sun  vector  in  the  spin  reference  frame  and  the  spin  rate  at  each  sample 
instant.  The  model  will  provide  a flag  indicating  whether  or  not  a pipper  command 
pulse  will  occur  during  the  next  sample  interval,  the  time  to  the  pipper  pulse  from 
the  current  sample  instant,  and  the  solar  aspect  angle  at  the  pipper  pulse  occurrence. 

DISCUSSION  OF  THE  MODEL 

A flowchart  of  the  preliminary  math  model  for  the  SSS  is  shown  in  Figure  E-l.  It 

accepts  from  the  main  program  the  spin  rate  and  the  sun  vector  in  the  spin  reference 

frame  at  the  sampling  instant.  To  simulate  the  asynchronous  pipper  pulse  within  a 

synchronous  simulation,  both  a flag  and  the  predicted  pipper  pulse  time  from  the 

current  sampling  instant  is  output  when  the  pipper  pulse  is  predicted  to  occur  within 

the  next  sample  period,  T . The  measured  solar  aspect  angle,  nominally  the 
‘ r samp 

angle  between  the  Z axis  and  the  sun  vector  and  specifically  the  angle  between  the  sun 
vector  and  the  sensor  Z axis,  will  be  computed  the  time  of  the  pipper  pulse, 
quantized,  and  output  to  the  main  program  along  with  the  flag  and  pipper  pulse  time. 

To  perform  these  computations,  the  model  must  be  initialized  with  the  angular  align- 
ment of  the  sensor  and  the  orientation  of  the  spin  vector  with  respect  to  the  vehicle. 

The  SSS  model  described  in  Figure  E-l  is  based  on  the  1/8/75  spin  sun  sensor  procure- 
ment specification  and  the;  assumptions  listed  in  Table  E-l.  Table  E-l  lists  the  initiali- 
zation data  required  by  the  model.  Table  E-2  lists  the  input  data  to  the  model  from 
the  main  program,  and  Table  E-3  lists  the  outputs  computed  by  the  model.  The  inter- 
mediate variables  internal  to  the  model  are  described  in  Table  E-4. 

The  notation  [0]j  denotes  a rotational  transformation  by  the  angle  6 about  the  j axis. 
Transformations  used  in  the  model  and  those  suggested  to  be  used  in  the  main  program 
to  obtain  SW,  the  sun  vector  in  the  spin  frame  at  the  sampling  instant  t,  are  defined 
in  Figures  E-2  through  E-7. 
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TABLE  E-l.  INITIALIZATION  DATA 


Variable 

Units 

Description 

aS 

Degrees 

Figure  E-6 

ss 

Degrees 

Figure  E-6 

Ys 

Degrees 

Figure  E-6 

“w 

Degrees 

Figure  E-7 

j eW 

Degrees 

Figure  E-7 

TABLE  E-2.  INPUT  DATA  FROM  MAIN  PROGRAM 


Variable 

Units 

Description 

1 SAMP 

Seconds 

SSS  Sample  Period 

ws 

Radians/ 

Second 

Spin  Rate 

v*> 

Sun  Vector  in  Spin  Reference  Frame  at 

Sample  Irfstant 

TABLE  E-3.  SSS  OUTPUT  TO  MAIN  PROGRAM 


Variable 

Units 

Description 

T . 

pipper 

Seconds 

Time  from  beginning  of  current  sample 
interval  to  time  of  pipper  pulse. 

SUN 

ANG 

Degrees 

Angle  between  Z sensor  axis  and  sun  vector 
(solar  aspect  angle). 

I . 

pipper 

Flag  set  to  1 if  proper  pulse  occurs  during 
current  interval.  Set  = 0 otherwise. 

TABLE  E-4.  INTERNAL  VARIABLES 


Variable 

Units 

Description 

iTv/sl 

Sensor  to  vehicle  transformation. 

Transpose  of  [Tg/y].  (See  Figure  E-6.) 

Vehicle  to  spin  transformation.  (See 

Figure  E-7. ) 

^Tw/s^ 

Sensor  to  spin  transformation 

fTW/S^  = ^TW/V^  ^s/v^ 

Boresight  vector  in  spin  frame.  (See 

Figure  E-5. ) 

Normal  to  sensor  plane  in  sensor  frame. 

(See  Figure  E-5.  ) 

8 

Radians 

Sensor  measured  value  of  complement  of 
the  sun  angle. 

Sensor  plane  normal  vector  in  spin  frame. 

ss 

Sun  vector  in  sensor  frame. 

ITER 

Iteration  counter. 

*S 

Unit  vector  along  sensor  Z axis. 

T P 

Seconds 

Value  of  Tpjpper  from  previous  iteration. 

8' 

Radians 

Value  of  B from  previous  iteration. 

The  method  of  solution  for  the  pipper  command  and  pulse  time  and  the  analysis 
performed  to  develop  the  model  are  documented  in  the  subsection  entitled  "Method 
of  Solution.  " The  analysis  showing  the  effect  of  the  quasi-stationary  spin  axis 
orientation  is  also  documented.  The  model  does  not  include  sensor  dynamics  or 
anomalies  in  the  optics. 

ASSUMPTIONS 

• The  orientation  of  the  spin  axis  to  the  vehicle  axes  is  defined  in  Figure  E-7 
and  is  constant.  Nominally,  the  spin  axis  and  the  Z vehicle  axis  are  aligned. 

• The  sensor  reference  frame  and  orientation  to  the  vehicle  axes  are  as  defined 
in  Figures  E-5  and  E-6.  Nominally,  the  Z sensor  axis  is  parallel  to  the  spin 
axis. 
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INERTIAL  FRAME 

X.  - in  celestial  equatorial 

plane  pointed  toward 
vernal  equinox 

Y.  - in  celestial  equatorial 

plane  at  90  degrees  right 
ascending 

Z.  - perpendicular  to 

equatorial  plane 
(nominally  toward  Polaris) 

TRANSFORMATION  FROM  INERTIAL 


ORBIT  FRAME 
Xg  - ascending  node 


Yg  - lies  in  orbit  plane 
and  90  degrees  from  Xg 

Zg  - perpendicular  to 
orbit  plane 
(right-hand  triad) 


TO  THE  ORBIT  ORIENTED  INERTIAL  FRAME  IS 


[v] ' ['.]  „ [”.]  z 


where  = right  ascension  of  ascending  node 


ij  = orbit  inclination 


Figure  E-2.  Definition  of  Transformation  from  the  Stellar  Inertial  Frame 
to  the  Orbit  Frame 
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J Y 
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This  is  a time  varying  transformation  driven  by  the  control  system  where 

Yq  = the  angle  from  Xp  about  ZQ  to  the  projection 
of  the  spin  vector  on  the  orbit  plane  XqYq 

= the  angle  between  the  spin  vector  and  the  orbit 
plane  normal 


Figure  E-3.  Transformation  from  the  Orbit  (O)  Frame 
to  the  Momentum  (H)  Frame 
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SENSOR  PLANE 
(e  = 01 


SENSOR  FRAME 


Y„  is  along  N,  the  normal  to 


= n plane  (p7ane  of 


the  narrow  field  of  view). 
-+ 

X$  is  along  B,  the  boresight 
which  is  the  intersection 
of  the  plane  e = 0 and  the 
plane  6=0. 

Zr  completes  the  right-hand 
orthogonal  triad. 


Figure  E-5.  Definition  of  GPS  SSS  Sensor  Reference  Frame 


I 
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M = [T  SPIN/VEHICLE j = 

This  is  a constant  transformation  where 

i-  = angle  between  the  Z„  axis  and  the  projection  of  the 
spin  vector  on  the  Xy  - Zy  plane 
= angle  between  the  spin  vector  and  the  Xy  - Zy  piane 


Figure  E-7.  Definition  of  GPS  Spin  Reference  Frame  and  Transformation 
from  Vehicle  to  Spin  Reference  Frame 


• The  sample  period  Tg^^p  will  be  0.05  sec  or  shorter  to  preclude  the  sun 
vector  from  entering  the  active  region  of  the  sensor  from  the  opposite  side 
of  the  vehicle  during  one  sample  period. 

• The  rate  of  angular  displacement  of  the  spin  vector  in  inertial  space  will  be 

less  than  0.01  deg/TSAMp  so  that  angular  displacement  of  th>-  :,pin  axis  during 
TSAMP  negligible.  This  will  allow  prediction  of  the  pipper  pulse  from 

the  spin  rate  and  geometry  without  accounting  for  spin  vector  displacement. 

• At  each  sample  instant,  the  simulation  main  program  will  provide  the  SSS 
model  with 

1.  Spin  rate  Wg  about  the  spin  axis  in  rad/sec,  and 

2.  the  sun  vector  in  the  spin  frame. 

• The  pipper  command  pulse  occurs  0.  25  deg/cos  B ahead  of  the  crossing  of 
the  sensor  (0  = 0)  plane.  The  repeatability  of  the  pulse  with  respect  to  9 = 0 
is  3 arc  minutes  for  repeated  crossings  at  constant  8. 

METHOD  OF  SOLUTION 

The  method  of  solution  is  to  predict,  from  the  sun  vector,  sensor  plane  unit  normal 
vector  and  the  spin  rate,  the  time  of  the  pipper  command  pulse,  and  the  measured 
solar  aspect  angle  at  that  instant.  It  is  assumed  that  the  spin  axis  is  quasi-stationary 
as  explained  in  the  next  subsection. 

Define  the  pipper  command  plane  with  unit  normal  Pg  to  contain  the  sensor  Z axis  and 
the  sun  vector  at  the  time  of  the  pipper  command  pulse.  The  unit  normal  Pg  can  be 
defined  in  terms  of  Ng,  the  normal  to  the  sensor  (9  = 0)  plane  shown  in  Figure  E-5. 

Ps  is  the  result  of  *he  rotation  of  Ng  0.  25  deg/cos  B about  the  Z sensor  axis,  where 
B is  the  complement  of  the  measured  solar  aspect  angle  also  shown  in  Figure  E-5. 

Ps  = (-0.  25/57.  3)*(1  / cos  8)  z Ng  (E-l) 

This  relationship  between  the  pipper  command  plane,  sensor  plane,  and  aspect  angle 
is  defined  in  Reference  E-l.  At  the  time  of  the  pipper  command  pulse,  the  sun  vector 
lies  in  the  pipper  command  plane;  therefore,  the  sun  vector  S and  the  normal  P are 
perpendicular  at  this  instant.  Thus, 


S • P = 0 


(E-2) 


The  solution  of  Equation  (E-2)  will  allow  the  prediction  of  the  pipper  pulse  time.  If 
sensor,  vehicle,  and  spin  reference  frames  all  have  Z axes  aligned,  the  solution  is  of 
closed  form.  However,  these  frames,  in  general,  will  not  be  aligned,  so  an  iterative 
solution  is  shown  in  Figure  E-2.  The  solution  of  Equation  (E-2)  will  be  done  in  the  spin 
reference  frame  so  that  only  rotations  about  one  axis  need  to  be  considered.  The  sun 
vector  in  the  spin  frame. 


SW~[1  W/lf*  f 1 II/O-I  ^ro/i^  si 


(E-3) 


where  the  transformations  from  the  inertial  frame  to  the  spin  frame  are  defined  in 
Figures  E-3,  E-4,  and  E-5,  will  be  input  from  the  main  program  at  the  sampling 
instant.  The  pipper  command  plane  unit  normal  in  the  spin  frame 

?W  = fTW/V]  ns/v]Tf(-0.25/57.3)  1/cos  6 Jz  Ng  (E-4) 


will  be  computed  in  the  model.  The  transformations  [T^^]  and  fT^y^]  are  defined 
in  Figures  E-7  and  E-8.  The  initial  value  of  Bis  computed  from  the  sun  vector  S 
at  the  sample  instant: 


W 


•■sin’1  <v  fTs/vmw/viT swi 


(E-5) 


At  the  instant  of  the  pipper  pulse,  T . , the  sun  vector  in  the  spin  frame  will  be 

1 pipper 

S„.(T  . ) = [W,,  T . j S,  , (E-6) 

W pipper  1 S pipper J Vv 

because  the  spin  frame  rotates  about  the  Z axis.  is  the  spin  rate  (about  the  Z spin 

axis)  and  T . is  measured  from  the  sample  instant.  Incorporating  (E-4)  and  (E-6) 
pipper 

into  (E-2), 


fW  T . ] S,„  • P,„  = 0 

S pipper  „ W W 
l Z, 


(E-7) 


Evaluating  Equation  (E-7),  one  obtains 


(sin  W0  T . )(S  P - S P ) + (cos  W_  T . )(S  P + S P ) = -S  P (E-8) 

S pipper  y x x y S pipper  x x y y 7.  7. 
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V (t+At)  Y (t) 

u U! 


P ( t+At ) Li  P(t) 


X (t+At) 


Figure  E-8.  Sun  Vector  Rotation 


which  has  the  form 


sin  9 B + cos  9 A = -K 


Equation  (E-8)  can  be  reduced  by  trigonometry  identities  to 

-S  P 

sin(W  T + 0)  = — - — ^ -a — 5 5 — 5-1 

S plpper  (S2P2  + s2p“  + s2p2  + s2p2)f 

XX  y y yx  xy 


where 


, S P +S  P 
- 1 xx  y y 
0 = sin  — - 

~\J (1  -S  2)(1  -P  2)  - 

v z z 


r (1  -S  2)(1  -P  2) 

z z 


T . = sin  ; — ; 

pipper  , r ^ 

h/(l-S  )(1-P  ) 


S P +S  P 
x x y y 

/ 2 2 
ln-s  )(i-p  ) 

z z 


Examining  the  second  term,  we  see  that  the  opposite  sign  is  needed  so  that,  for  positive 

spin  rate  a nd  P , P , S , S positive,  the  time  T . will  be  positive  for  the  time 
x y x y r pipper 

required  for  the  sensor  to  rotate  so  that  the  sun  vector  is  perpendicular  to  the  plane 

normal  P (Figure  E-9).  Examining  the  first  term,  we  see  also  that  the  sign  is  wrong 

because  the  sensor  must  be  rotated  farther  to  cross  the  sun  vector  when  S and  P are 

z z 

of  the  same  sign  than  when  they  are  either  zero  or  of  opposite  signs  (Figure  E-9).  Thus, 
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COMMAND 

PLANE 


Because  the  sensor  frame  is  not  aligned  to  the  spin  frame,  the  measured  8 is  recomputed 
for  the  sun  vector  in  the  sensor  frame  at  T . Equation  (E-10)  is  then  iterated  with 

pp*--  1 

Phased  on  the  new  value  of  B.  This  iterative  process  can  be  repeated  until  the  Tp.pper 

and/or  B converge  to  acceptable  tolerances  or  for  a maximum  number  of  iterations. 

if  -r  > T , the  pipper  command  will  occur  following  the  next  sampling.  If 

pipper  samp 

T <0,  the  pipper  command  has  already  occurred  in  a previous  sample  period, 

pipper 


There  is  also  a random  error  which  must  be 
repeatability  of  the  pipper  command  at  6=0. 
same  for  repeated  passes  through  any  angle 


added  to  T . It  is  the  3 arc  minute 

pipper 

This  repeatability  is  assumed  to  be  the 
B . The  random  time  error  will  be  modeled 


by  a nominally  distributed  random  number  representing  the  uncertainty  in  the  time  of 
the  pipper  pulse.  The  1 sigma  value  of  this  time  uncertainty  is 


a 


rand 


1/20  _l_ 

57.3  W 


sec 


If  the  spin  rate  is  variable,  a random  number  generator  which  has  a variable  standard 
deviation  will  be  needed. 
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SIMULATION  ACCURACY  WITH  ASSUMPTION  OK 
QUASI- STATIONARY  SPIN  AXIS 


The  GPS  SSS  is  designed  to  operate  in  the  spin  rate  range  of  5 to  100  rpm.  To  model 
the  sun  crossing  pipper  pulse  time  hut  maintain  a relatively  long  sample  period  T^ 
it  is  necessary  to  predict  the  time  within  the  sample  period  during  which  the  pipper 
pulse  occurs.  It  is  assumed  that  the  rate  of  change  of  the  spin  vector  orientation  is 
sufficiently  small  (loss  than  .01  deg/T^  p)  so  that  the  angular  motion  of  the  sensor 
with  respect  to  the  sun  vector  can  be  approximated  by  the  angular  motion  about  the  spin 
axis  location  at  the  last  sample  instant.  Then  the  prediction  of  pipper  pulse  time  is  a 
straightforward  computation  from  the  spin  rate  and  the  sun  vector  in  the  spin  frame  at 
the  last  sample  instant.  If  the  spin  vector  orientation  will  be  changing  faster  than  .01 
deg/T^^p,  then  the  model  needs  lobe  changed  to  use  an  iterative  method  or  to  extra- 
polate the  sun  vector  to  the  spin  frame  at  the  instant  of  pipper  crossing. 

The  maximum  spin  axis  orientation  rate  according  to  the  Reference  E-3  precession 
control  mode  is  . 186  deg  per  jet  firing.  The  duration  of  firing  is  .111  sec  and  firing 
occurs  once  per  revolution.  Since  the  precession  equations  are  derived  from 

::t  = + w x u 

or 

T = W \ hk 

T = W h 

* y 

T = -W  h 

y x 

the  jet  torque  causes  rate  rather  than  acceleration.  The  precession  rate  will  be 
.186  deg/. Ill  see  =1.66  deg/sec. 

This  rate  of  change  of  the  spin  axis  orientation  does  not  occur  as  a continuous  function, 
but  the  rate  is  present  for  about  60  deg  out  of  each  rotation.  With  a sampling  period 
of  .005  sec,  the  error  in  spin  axis  orientation  knowledge  at  the  predicted  pipper  pulse 
lime  would  be  about  . 008  deg  in  the  worst  case.  Worst  case  will  occur  when  the  pipper 
pulse  occurs  at  the  end  of  the  sample  interval.  This  will  be  an  error  of  .008  deg  in 
the  aspect  angle  measurement  with  negligible  error  in  the  pipper  pulse  time. 
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APPENDIX  F 


REACTION  JET  MODEL  FOR  GPS  INDEPENDENT  ANALYSIS 


This  appendix  presents  the  model  for  the  GPS  reaction  jets.  This  model  is  based  on 
the  model  used  in  the  3XATS  simulation.  Both  ATS  and  GPS  use  Rocket  Research 
Corporation  monopropellant  hydrazine  thrusters.  Any  data  that  was  not  available  from 
Rockwell  was  extracted  from  ATS  documents. 

MODEL  DESCRIPTION 

The  model  consists  of  an  initialization  section  and  a normal  operation  section.  The 
initialization  section  performs  one-time  computations  and  clears  outputs  and  counters. 
The  normal  operation  section  receives  jet  commands  from  the  controller  and  outputs 
total  jet  force  and  total  jet  torque  on  the  spacecraft.  The  following  are  system  charac- 
teristics included  in  the  model: 

• Thrust  increases  with  catalyst  temperature  which  increases  when  the  jet  is 
firing  and  decreases  when  the  jet  is  not  firing.  The  time  constant  for  heating 
is  less  than  that  for  cooling. 

• No  thrust  is  developed  if  the  jet  command  is  less  than  the  minimum  on  time. 

• Thrust  start  and  stop  times  lag  the  jet  command  start  and  stop  times  due  to 
valve  response  time  and  hydrazine  burn  response. 

• Eighteen  thrusters  are  used. 

• Thruster  torque  depends  on  variable  spacecraft  eg  location. 

The  following  are  jet  system  characteristics  not  included  in  the  model: 

• Dynamics  of  the  thrust  rise  at  turn-on  and  thrust  decay  at  turn-off  must  be 
considered  because  of  the  impact  on  computer  run  time.  The  static  model 
will  be  adequate  for  all  but  difficult  limit  cycle  cases.  These  limit  cycle 
cases  should  be  studied  using  a single-axis  simulation  and  correlated  with 
static  model  results.  A dynamic  model  can  be  provided  if  necessary. 


• Changes  in  thrust  due  to  supply  pressure  reduction  with  fuel  consumption  is 
not  modeled.  The  supply  pressure  changes  will  be  very  slow.  Stability  should 
be  analyzed  with  both  high  and  low  supply  pressure  by  changing  the  thrust  level 
initialization  inputs. 

• Dynamic  changes  in  thruster  position  and  angular  alignment  are  not  considered. 
The  spacecraft  is  very  stiff.  Cases  of  extreme  position  and  angular  alignment 
should  be  run. 

The  jet  model  flow  diagram  is  shown  in  Figure  F-l.  The  left  arrow  (— ) within  boxes 
indicates  replacement,  i.  e. , the  expression  on  the  right  replaces  the  variable  on  the 
left.  Mathematical  subscripts  are  used.  Symbols  used  in  the  flow  diagram  are  defined 
in  the  sections  below.  Comments  appear  outside  the  boxes  to  describe  the  operations 
within  the  boxes. 

INITIALIZATION  DATA 

Table  F-l  lists  the  initialization  parameters,  nominal  values,  tolerance,  and  source 
of  data.  This  data  is  to  be  input  at  the  start  of  each  run. 

INPUT  AND  OUTPUT  DATA 

Table  F-2  lists  I he  input  data  that  the  jet  model  expects  from  the  calling  program  each 
time  the  jet  subroutine  is  executed  in  the  normal  mode. 

Table  F-3  lists  the  output  data  from  the  jet  model  subroutine.  Other  subroutine 
variables  may  be  outputs  if  the  user  desires. 

INTERNAL  VARIAHLES 

Internal  variables  used  in  the  flow  diagram  are  defined  in  Table  F-4. 
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( ENTER  > 
^SUBROUTINE  JETSj 


NOTE  LEFT  ARROW  (• » INDICATES 
REPLACEMENT  OPERATION 


Figure  F-l.  Flow  Diagram  of  Reaction  .Jet  Model 


INCREMENT  TIME  COUNTER 


CLEAR  TORQUE  FORCE  AND 
JET  COUNTER  OUTPUTS  FOR 
UPDATE 


UPDATE  HISTORY 
OF  JET  COMMANDS 


PUT  LATEST  JET 
COMMAND  INTO  HISTORY 


(JET 

TURNING  ONI 


(JET  COMMAND 
REMAINING  OFF 


i JF  T commano 

REMAINING  ON  ) 


(JET  TURNING  OFF  I 


TXJET 

0 

TYJET 

- 

0 

TZJET 

0 

FXTOT 

0 

FYTOT 

- 

0 

FZTOT 

- 

0 

JETCT 

0 

p 


> 0 


(JET  WILL 

:_0  STAY  ON  FOR  f 

PART  OF  INTERVAL.) 

(JET  IS 
OFF.) 


DECREASE  JET  THRUST 
POTENTIAL  DUE  TO 
CATALYST  COOLING. 
SET  JET  THRUST  TO  0. 


FSAV,  • FSAV,  * (FMIN;  FSAV,)  * COFF, 
FJET,  - 0 


Figure  l’-l.  Flow  Diagram  of  Reaction  Jet  Model  (continued) 


I ABI..E  1-1.  INITIALIZATION  DATA  FOB  JET  MODEL 


Parameter  Subscript  Units 


Nozzle  X 1 

coordinate  2 


Nominal 


2.  833 
2.  833 
2.  833 
2.  833 
-2.833 
-2. 833 
-2. 833 
-2.  833 
2.  667 
2.  667 
2. 667 


Source 


Tolerance 


Source 


Nozzle  Y 
coordinate 


Not 

available 


Scaled  from  Rockwell  drawing  v 505-900003,  Rev.  C,  1+27-75. 

Rockwell  GPS  PDR  handout. 

Coordination  meeting  notes,  3/20/75. 

Rockwell  Specification  MC  321-0009,  10/24/74. 

Rockwell  Specification  MC  321-0011,  10/24/74. 

GPS  0.1  I.llE  Rl'A,  Qualification  Test  Report,  CDRI,  Hem  A005,  March  31, 
1976. 

GPS  5.0  l.BF  RKA,  Qualification  Test  Report,  CDRI,  Item  A005,  March  31. 
1076. 

Pairs  match  to  _+  5%,  3c.  Pairs  are  9 and  14,  10  and  13,  12  and  15,  11  and 
16,  3 and  6,  4 and  5,  1 and  7,  2 and  8,  6 and  7,  2 and  4,  1 and  3,  5 and  8, 

9 and  15,  10  and  16,  12  and  14,  11  and  13. 


T ABLE  F-l.  INITIALIZATION  DATA  FOR  JET  MODEL  (continued) 


Flow 

Definition 

Diagram 

Symbol 

Dimension 

Type 

Units 

Total  number  of  jet 
operations 

JETCT 

Integer 

Number  of  operations  for 
each  jet 

JETCTR 

1 8 

Integer 

Accumulated  on  time 
for  each  jet 

TOTJT 

18 

Real 

Seconds 

X component  of  total  jet 
force  on  spacecraft 

FX  TOT 

Real 

Pounds 

V component  of  total  jet 
force  on  spacecraft 

FYTOT 

| 

Real 

Pounds 

Z component  of  total  jet 
force  on  spacecraft 

FZ  TOT 

Real 

Pounds 

Total  jet  torque  about 

X axis 

TXJET 

Real 

Foot 

pounds 

Total  jet  torque  about 

V axis 

TV  JET 

Real 

Foot 

pounds 

Total  jet  torque  about 

Z axis 

TZ  TFT 

Real 

Foot 

pounds 

* 


TABI.E  F-4.  INTERNAL  VARIABLES  FOR  JET  MODEL 


Flow 

Definition 

Diagram 

Symbol 

Dimension 

Type 

Units 

Thrust  direction  cosine 

CACB 

18 

Real 

... 

Thrust  direction  cosine 

CASH 

18 

Real 

... 

Thrust  direction  cosine 

SA 

18 

Real 

... 

Thrust  rise  increment 

CON 

18 

Real 

... 

Thrust  decay  increment 

COFF 

18 

Real 

... 

Potential  jet  thrust 
from  previous  cycle 

FSAV 

18 

Real 

Pounds 

Predicted  time  for  jet 
thrust  to  go  off 

PTOFF 

18 

Real 

Seconds 

Quantized  on  delay 

M 

18 

Integer 

Cycles 

History  of  jet  commands 

JET2 

18x20 

Integer 

... 

Time  since  start  of  run 

TY  M 

Double 

Seconds 

precision 

Quantized  minimum  on 
time 

ICMIN 

Integer 

Cycles 

Incremental  jet  on  time 

DON 

Real 

Seconds 

Actual  jet  thrust 

FJET 

in 

Real 

Pounds 

X component  of  jet  force 

FX 

18 

Real 

Pounds 

Y component  of  jet  force 

KY 

18 

Real 

Pounds 

Z component  of  jet  force 

FZ 

18 

Real 

Pounds 

X component  of  nozzle 
moment  arm 

XI. 

18 

Real 

Feet 

V component  of  nozzle 
moment  arm 

YL 

18 

Real 

Feet 

Z component  of  nozzle 
moment  arm 

ZL 

18 

Real 

Feet 

X component  of  jet  torque 

TXJ 

18 

Real 

Foot  pounds 

Y component  of  jet  torque 

TY  J 

18 

Real 

Foot  pounds 

Z component  of  jet  torque 

TZJ 

18 

Real 

Foot  pounds 

Jet  index 

i 

Integer 

General  index 

.1 

Integer 

t 
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APPENDIX  G 


GPS  DEPLOYMENT  MODEL 


The  following  equations  will  be  used  for  simulation  of  the  GPS  rotational  dynamics 
during  deployment.  These  equations  will  be  incorporated  in  the  one-body  simulation 
to  investigate  the  effect  of  deployment  on  the  pitch  and  roll  jet  channels. 

The  following  primary  assumptions  were  made  in  developing  these  equations: 

1.  Panels  are  treated  as  point  masses. 

2.  Mass  of  the  boom  is  neglected. 

3.  A "linear"  deployment  linkage  is  assumed. 

R OT  ATIO  N A L 1 )Y  N A M I C.S 

Referring  to  Figure  G-l  and  Appendix  O for  notation,  the  angular  momentum  during 
deployment  is  given  by 


II  1 u.>  + Lm.(r  f r.  + p.)  x p. 
0 0.10  1 1 1 
1 


where 


I J - m M ( r ,i  ) + r m.MIr.  + p.,  r.  + p.) 
o o 0 0,11111 

l 


The  rotational  differential  equations  of  motion  are  given  by 


1 - I 'i,  + I (i)  4 ;m.(r  x p.)  + fm.(r  + r.  + p.)  x p.  + id  x H 

0000  .10  1 .10  1 1 l o 

i i 


The  deployment  linkage.  Figure  G-2,  is  such  that  6 - 2a . With  this  linear  relationship, 
derivatives  appearing  in  Equations  (G-l)  through  (G-3)  can  be  evaluated  as  follows: 


p.  < (-1 ) (a  + b)  sin  <i. 

’ I 1 

I (!'  - a)  cos  a . 


t 
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APPENDIX  H 

SOLAR  PRESSURE  TORQUE  AND  FORCE  MODEL 

This  appendix  summarizes  the  preliminary  solar  pressure  force  and  torque  model 
developed  for  the  digital  simulation  evaluation  of  GPS  stability  and  control  analysis. 

The  solar  force  results  fro?n  the  photo  momentum  transfer  to  the  vehicle  surfaces 
exposed  to  solar  illumination.  The  magnitude  and  direction  of  the  solar  force  developed 
on  a given  vehicle  surface  are  functions  of  the  vehicle-sun  attitude,  the  geometry  of 
the  surface,  and  the  physical  characteristics  of  the  surface.  The  solar  torque  will  be 
developed  if  the  vehicle  center  of  mass  is  not  coincident  with  the  center  of  total  net 
force. 

The  model  described  in  this  appendix  was  based  on  the  vehicle  geometrical  configur- 
ations deduced  from  the  Rockwell  drawing  No.  V505-900003  titled  "Spacecraft  General 
Arrangement  Baseline  Configuration  - GPS,  " January  27,  1975. 

VEHICLE  GEOMETRY  FOR  SOLAR  TORQUE  EVALUATION 

The  GPS  spacecraft  configuration  has  been  approximated  by  surfaces  of  simple  geo- 
metry to  facilitate  the  evaluation  of  solar  force  and  torque.  This  simplified  geometry 
consists  of  five  flat  surfaces,  six  cylindrical  surfaces,  one  segment  of  a cone,  and 
fourteen  complete  cylindrical  surfaces  for  the  representation  of  vehicle  main  body, 
solar  panels,  thruster  cone,  and  antenna  in  on-orbit,  configuration.  The  stowed  vehicle 
during  the  spinning  phase  is  approximated  by  a complete  cylindrical  surface.  The 
geometries  of  the  simplified  vehicle  configurations  are  shown  in  Figure  11-1. 

Total  solar  pressure  force  and  torque  experienced  by  the  vehicle  can  be  evaluated  as 
the  vector  sum  of  the  forces  and  torques  exerted  on  the  approximating  surfaces. 
Nominal  geometries  including  size,  location,  and  orientation  of  the  surfaces  for  the 
GPS  vehicle  are  summarized  in  a later  section  of  this  appendix.  These  data  are  based 
on  the  Rockwell  drawing  No.  V505-900003  titled  "Spacecraft  General  Arrangement 
Baseline  Configuration  - GPS,"  January  27,  1975.  The  computations  involved  in  the 
evaluation  of  the  force  and  torque  contributed  by  these  elementary  surfaces  are 
described  in  the  sections  that  follow. 


Solar  Pressure  on  Flat  Surface 


A flat  surface,  defined  in  terms  of  the  surface  normal  N and  geometrical  center  C, 
develops  a solar  pressure  force  F and  torque  jr  about  vehicle  center  of  mass  as  given 
in  the  following  equations: 

F = -F^N  + Ftf(N  • S)N  - S]/  !(N  • S)N  - S 1 (H-l) 

l=CxF  (H-2) 

where 

S = LOS  vector  to  sun 

F = [(1  + so)  cos"  -j  + rr  p ( 1 - s)  cos  \|i]  PA 
n i 

F = f(l  - sp)  cos  il  sin  ijr  ] PA 

-7  9 

P = solar  pressure  = 1x10  lb/ft 

P = fraction  of  incident  radiation  reflected 

s = fraction  of  reflected  radiation  that  is  specular 

C = position  vector  of  surface  geometrical  center  from  vehicle  center  of  mass 


cos  - (N_  • S) 

Note  that,  for  the  case  N = S, 

| (N  • S)N  - S | = 0 

This  portion  of  computation  in  Equation  (11-1 ) should  be  bypassed.  This  is  when  there 
is  no  tangential  force. 
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From  the  geometry  of  Figure  11-2,  the  location,  orientation  and  size  of  a cylindrical 
surface  can  be  completely  specified  in  terms  of  the  following  parameters: 

C_  = position  vector  of  the  geometrical  center  of  the  surface 
N = normal  vector  at  C (orientation) 

= tangent  vector  at  C (orientation) 
r = radius  of  the  cylinder  (size) 

9q  = angular  width  of  the  cylinder  (size) 

H = length  of  the  cylinder  (size) 

The  evaluation  of  solar  force  on  a cylindrical  surface  can  be  carried  out  by  the  sub- 
division into  n flat  surface  elements  as  shown  in  Figure  H-3.  The  i1*1  flat  surface 
element  is  defined  in  terms  of  the  location  of  its  geometrical  center  C.  and  the  normal 
vector  N.  which  are  related  to  the  cylinder  parameters  as 


Figure  11-2.  Surface  Numbering  Notation 


C.  = C - rN  + rN. 

N = cos  @.N  + sin  P.T 

— 1 l—  l— 


where 


9i  = '90/2  " A/2  +iA*  1 = 


0 In 


The  area  of  i1^1  surface  element  is 


A.  = r • A • i 

l 


(H-3) 

(H-4) 

(11-5) 


(H-6) 


The  solar  forces  and  torques  developed  on  these  flat  surface  elements  can  be  computed 
using  Equations  (11-1)  and  (11-2).  Let  the  forces  and  torques  be  denoted  as  F.  and  t.. 

The  total  force  and  torque  developed  by  the  cylindrical  surface  are  Fc  and  evaluated  as 


F = I F. 

-c  i=i-> 


T = T T 

=l~C 


(H-7) 


(11-8) 


For  the  case  of  a full  cylindrical  surface,  such  as  the  antenna,  the  cylinder  is  defined 
in  terms  of  following  parameters: 

C = position  vec'or  of  the  center  of  the  cylinder  from  vehicle  center  of  mass 
L = unit  vector  along  the  cylinder  axis 
r = radius  of  cylinder 
1 - length  of  cylinder 

The  solar  force  F and  torque  t about  vehicle  center  of  mass  become 


(11-9) 

(H-10) 


1 


F = [(1  - So)  sin  0 cos  e]PA 
A = projected  cylinder  area  = 2£r 


cos  e = (S  • L) 


D = C + rN 


Again,  the  computations  associated  with  S - (S  • L)Lj  should  be  bypassed  for  S = L. 
This  is  when  theme  is  no  transverse  force. 


Solar  Force  on  Conic  . urfac e 

From  the  geometry  of  Figure  11-4,  the'  location,  orientation,  and  size  of  a conic  surface 
can  be  completely  specified  in  terms  of  the  following  parameters: 

N = unit  vector  along  cone  axis 
— c 

r = radius  of  circle  on  top 

r,  = radius  of  circle  on  bottom 
b 

o = center  of  circle  on  top 
—a 

h = height  of  conical  surface 

The*  evaluation  of  the  solar  force  on  a conical  surface  can  be  carried  out  through 
approximations  using  flat  surface  elements.  Each  surface  element  is  a trapezoid  which 
can  be  represented  by  a surface  normal  N.  and  the  location  of  its  geometrical  center  C. 
as  depicted  in  Figure  11-4.  These  parameters  art'  related  to  the  defining  parameters 
of  the  conical  surface'  as 


N.  = cos  0 cos  n.  X + cos  0 sin  n.  Y + sin  o N (H-l  1 ) 

— l c l — c c i — c c— c 


C.  = A.  + (2r  + r )/ (3rb+  3r  ) f B.  - A.] 

— l—l  b 1 1 —i  — iJ 


(H-l  2) 


where  X(  and  arc  unit  vectors  that  make  (X(  - Yf  - N ) a right-handed  orthogonal 


triad.  For  convenience,  let 


Y = (N  x Xj/|N  x X I.  X = - Y = -X_ 
— c -c  —s"  -c  —s'  SBC. 


a = i6  - 6/2,  i = 1,  2, 


6 = rr/n 


88 


GEOMETRICAL  CENTER  OF 
TRAPEZOID: 


tel 

I2  rb  + ra  \ 
\3rb  + 3ra  / 


Figure  11-4.  Geometry  of  Thruster  Cone 


l.  = O + r cos  r,  .X  + r sin  nr.  Y 
■l  — n a l—  c a l— c 

I.  = (O  - t,  N ) + r,  cos  a . X + r,  sin  nr  "S 

•l  —a  1— c b l— c b c- 

c = tan'1  (<rb  - r^/h  ^ 

1 * ra  + 8 (rb  ' ra)A/rb  ' ra>2  + h2 

, = distance  between  A.  and  B. 

1 —l  —l 

' distance  between  AL  and  TV 


The  area  of  the  trapezoid  is 


A.  = 4 (r.  + r,  ) * (H-13 

l l 1 b 1 

The  solar  force  and  torque  contributed  by  each  of  the  trapezoidal  surface  elements  can 
then  be  computed  using  Equations  (I  I - 1 ) and  (H-2)  with  the  total  force  and  torque  contrib 
uted  by  the  thruster  cone  evaluated  in  the  same  manner  as  Equations  (H-7)  and  (H-8). 


Vehicle  Nominal  Attitude  and  Shading  Considerations 

The  GPS  spacecraft  on-orbit  attitude  is  sun  and  earth  oriented.  The  pitch  and  roll 
control  allows  the  vehicle  yaw  axis  to  be  aligned  with  local  vertical.  The  vehicle  yaw 
motion  is  such  that  the  solar  array  surface  normal  can  be  driven  to  be  aligned  with  sun 
lino  for  maximum  solar  exposure.  Attitude  accuracy  requirements  are  0.  5 deg  for 
roll  and  pitch  and  3 deg  for  yaw. 


In  order  for  the  vehicle  main  body  to  cast  shadow  on  the  solar  array,  the  minimum  yaw 

attitude  error  must  be  * . 

s 

*3  = tan  **  tan  1(^|)  «=  27.  5 deg 

where 

f = clearance  between  main  body  and  solar  array 
b = length  of  main  body 

Similarly,  the  main  body  will  not  be  shaded  by  the  solar  array  panels.  Also,  the  solar 
array  attitude  is  driven  to  within  5 deg  of  maximum  exposure,  lienee  it  can  be  assumed 
that  the  solar  array  is  completely  free  of  shadow.  Since  the  sun  vector  can  be  any- 
where within  the  pitch  plane,  the  shading  of  the  thruster  cone  and  the  flat  surf  ce 
interfacing  the  thruster  cone  need  to  be  considered.  The  shading  of  the  antenna  is 
ignored  due  to  the  relatively  small  areas  involved. 

Shading  of  Thruster  Cone  by  Surface  P-,  --The  shading  of  the  thruster  cone  can  be 
considered  in  terms  of  the  i1^1  flat  surface  element  described  earlier.  Referring  to 
the  geometry  of  Figure  11-5,  the  projection  of  point  H.  on  the  surface  P^  is  1)^ . Any 
point  located  on  the  surface  P satisfies  the  vector  equation 

<j?i  - sw  -P1  = 0 (n-i4> 

Also,  the  vector  (11.  - B.)  must  be  pointed  along  the  sun  line',  i.e., 

— l—i 


(I).  - B.)  = rS 

— l —l  — 


(11-15) 


where  r distance  between  B.  and  I)..  The  proioclion  of  vector  (B.  - A.)  on  P.  along 

—l—i  J — l — l 1 b 

S is  the  vector  (1).  - A.)  which  intersects  the  edge  of  P.  at  E..  lienee  E.  satisfies 

— —i—i  1 — i —i 

the  vector  equations 


(E. 

— l 


%-iW 


2 

r 

c 


(H-16) 


(E.  - A.)  = t(l).  - A.) 

—i  —l  —i  —l 


(11-17) 


where  t - ratio  of  distance  between  E.  and  A to  that  between  I)  and  A.. 

— i — i —l  —l 


!)1 


I 


C„1  (3)  ■ H 


Xc  = (0,  1,  0) 


jag 


Figure  II— 5.  Shading  of  Thruster  Cone  by  Plano  Pj  Illustrated 
in  Terms  of  Elementary  Surface  A.  - B. 


Solution  of  the  vector  equations  for  D.  and  E.  is  included  in  this  appendix.  Once  the 
point  E.  is  determined,  the  shade  of  E.  on  the  surface  element,  denoted  as  the  point  A', 

—l  —l  r — j 

can  be  computed,  i.  o.  , 


A!  = A.  + U /i  )(B.  - A.) 
— i —l  AE  AO  —l  —l 


where 


X „ , = distance  between  A.  and  E. 

AE  —l  —l 

X . _ = distance  between  A.  and  D. 

AD  — i — i 


To  evaluate  the  solar  force  contributed  by  the  unshaded  portion  of  the  i™  surface  ele- 
ment, Equations  (11-11),  (11-12),  and  (11-13)  can  be  used  with  A.  replaced  by  A!.  It  may 
happen  that  the  entire  i1*1  surface  element  is  in  the  shadow  of  P . The  condition  for 

this  situation  is  that  the  distance  between  C’  and  D.  must  be  less  than  the  radius  of 

— PI  —i 


surface  P;'s  circular  edge  r , i.e.. 


<Cpl-p.).(Cpl-p.)<rc 


(H-19) 


If  the  inequality.  Equation  (H-19),  is  valid,  the  computations  associated  with  the  itn 
surface  element  can  be  bypassed  since  it  contributes  nothing  to  the  total  force  and 
torque.  Also,  it  may  happen  that  the  entire  i1*1  surface  clement  is  not  shaded  by  P , 

This  is  the  case  when  the  sun  vector  has  positive  projection  along  the  normal  to  P , i.  e. , 


cos  9 > (S  • N ) > 0 

c — —PI 


(H-20) 


If  the  inequality.  Equation  (11-20),  is  valid,  the  point  A!  is  set  to  be  identical  to  A.. 

Shading  of  Surface  P^  by  'Thruster  Cone- -Since  the  circle  defining  the  open  end  of  the 

thruster  cone  is  parallel  to  the  surface  P , the  projection  of  this  circle  on  P^  is 

always  a circle  of  identical  size  as  shown  by  the  parallelograms  of  Figure  H-6.  The 

geometries  of  the  shadow  on  surface  Pj  for  four  different  sun  line  configurations  are 

considered  for  the  following  discussions.  The  shadow  of  thruster  cone  on  P^  will  be 

characterized  as  a flat  surface  of  area  A^  with  normal  N and  geometrical  center  C^. 

Equations  (11-1)  and  (11-2)  are  then  used  to  compute  the  contribution  of  solar  force 

F and  torque  t as  if  A were  not  shaded.  Hence,  the  effective  solar  force  and  torque 
— s — s s 

contributed  by  the  shaded  surface  P is 

F = F (0)  - !•' 

-PI  -PI  -s 


r - t - t 

-PI  -PI  -s 


where 


F ^ and  t p^  are  force  and  torque  computed  for  a completely  unshaded 


surface  P . 


Let  the  center  of  the  circle  where  thruster  cone  joins  P^  be  C.  , and  let  the  center  of 

the  circular  shadow  of  the  end  of  thruster  cone  on  P„  be  (1  . Position  of  Cl  can  be 

1 — b — b 

solved  from  the  vector  equation 


- ca>  • np1  o 


C.  - (C  f h N . ) = -X 
— b —a  —PI  s 


lere  X = length  of  [C^  - (C  ^ + h Np^)l. 
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Case  2: 


1 2 

A,  = 77  r |2g  - sin  2 6] 

12a 

<£-,  " £V  2/3  r.l  sin'  p 

— 1 — a i . B- cos  B sin  6 

ab 


1 2 
At,  = 77  (‘2r  sin  R + 2 r,  sin  R)  a . sin“B 

n on 


C = fC  + l.(r  cos  B)]  - L h (2r  + r )/(3r  + 3r  ) 

— 2 —a  — a — b a b a 

= C + L \(r  cos  B)  - i , sin^B  (2r  + r )/(3r.  + 3r  )] 
—a  — a ab  b a b a 


1 2 ro  .2 
A.,  =r  r [2B  - sin  B 
— 3 2 b 


C = C f L tr  sin  B)/  (B  - cos  8 sin  ft) 
— 3 — b — 3 1) 


Total  area  of  shadow: 


As5Al+A2+2nVA3 


Centroid  position  of  A : 

s 


S.  f IAiS1  + A2£2‘5’"'b2cb-A.!c3] 

s 
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Centroid  of  A : 

s 

Cs=frAA  +A5£5+  *<£61 

s 

where 

cos  28c  =[(Ej  - Ca)  • (E2  - CJ/[(Ej  - CM  • (Bj  - Cg)] 

E.  = F.  +\/r  2 - r2  (G.  - F.)/  I (G.  - F.)l  i = 1,2 

— i — i V c a —l  —i  —i  —l 

F.  = C + r cos  P L+(-l)‘+1  r sin  BM 

— i — a a — a — 

G.  = C,  + r,  cos  8 E + (-l)i+*  r,  sin  8 M 

— i — b b — b — 

M = NP1  k 

SIMPLIFIED  GPS  VEHICLE  GEOMETRY  FOR 
SOLAR  FORCE  AND  'I'ORQUE  EVALUATION 

Planar  Surfaces 

The  five  planar  surfaces  include  P where  thruster  cone  is  located,  P,?  where  antennae 

are  located,  P and  P adjacent  to  the  solar  panels  in  the  directions  of  +Y  and  -Y_ , 

.»  4 — I > — J j 

respectively,  and  the  plane  introduced  at  the  open  end  of  the  thruster  cone  to  simplify 

the  complex  geometry  involved.  The  vehicle  body  frame  (B-frame),  as  defined  in  the 

Rockwell  drawing  for  spacecraft  general  arrangement  configuration,  has  its  Y^;>  axis 

parallel  to  the  line  of  the  solar  array  drive,  the  axis  along  the  AKM  thruster 

centerline,  and  the  X axis  to  complete  a right-handed  orthogonal  triad.  The  origin 

of  the  B-frame  is  located  at  the  center  of  the  open  end  of  the  thruster  cone  (not  the 

AKM  nozzle)  as  shown  in  Figure  1 1 - 7 . Let  each  plane  be  defined  in  terms  of  its  surface 

normal  N.  location  of  surface  centroid  C.,  and  area  of  surface  A.  in  B-frame  coordin- 
— i —l  i 

ates. 


Plane  P : 


Plane  Pn: 


(Sp2>  = 


0 

0 

61.  38 


P2)  1 ° 


AP2  = AP1 


Plane  P^ : 


^P3>  = 


0 

21.  8 

42.  63 


(N  p3)  - 1 


Ap3  = 37.  5 x 41. 8 
= 1567. 5 


Plane  P : 

4 


^P4>  " 


0 

-21.  8 

42.  63 


<NP4)-  -1 


A„.  = 37,  !i  x 41.8 
P4 

= 1567. 5 


Plane  P,: 
3 


(CP5>  - 


<*P5> 


Ap5  = n (14.  6) 
= 669.7 


Cylindrical  Surfaces 

The  six  segments  of  cylindrical  surfaces  include  two  main  body  surfaces,  and  C9, 
in  the  directions  of  Xp  and  -Xp,  respectively,  and  two  solar  panel  segments  for  each 
of  the  two  solar  arrays.  Each  segment  of  cylindrical  surface  is  defined  in  terms  of 
the  location  of  the  geometrical  center  C,  the  normal  and  tangent  vectors  at  C,  N,  and 
T,  respectively,  the  radius,  length,  and  (he  angular  width  of  the  cylinder  segment, 
r,  l,  and  9,  respectively.  The  cylindrical  surfaces  are  defined  as  follows: 


Cylindrical  Surface  C : 


li 

.80.2 

(£ci)s 

0 

42.  63 

r 

rci  = 32 

°-  xci 

(Nd* 


(Ici>B  1 


Ji 


Cylindrical  Surface  : 
_ B 

&C2'  = 


-30.2 

B 

'-1 

B 

" 0 

0 

<NC2>  = 

0 

(T  ) = 

- C2' 

-1 

_ 42.63 

. 0 . 

0 _ 

rc2  = 32.0,  iC9  = 37.5,  0C2  = 86.4  deg 

Cylindrical  surfaces  C^  through  Cg  are  solar  panel  surfaces.  In  the  following,  the 
normal  and  tangent  vectors,  N and  T,  will  be  defined  in  terms  of  the  panel  fixed  frames. 
Their  B-frame  definition  will  be  obtained  through  a transformation  as 


B 


B 


'iWc'  “"‘"V'Wic1 


The  positions  of  centroid  of  the  cylindrical  surfaces  will  be  defined  in  P-frame  with 
respect  to  origin  of  P-frame  as 


t(cc)  - <ocn 


Their  P-frame  definitions  are  then  obtained  as 

B P P 


B 


(Cc)  = Tnp[(Cc)  - (Oc)]  - (Oc) 


Each  cylindrical  segment  is  individually  defined  so  that  misalignments  between  segments 
of  arrays  and  between  arrays  can  be  simulated.  P-j:  PANEL  FIXED  FRAME 


Nominally 


<03)  = (0,  21,9,  42.  38)'1 


0 = 98.  6 deg 
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Nominally , 


— P5  = -T5 


% 


defines  Xp5  wilh  respect  to  XR 

S5  = rC5  ' rC5  C°S  <0C5/2) 
dr  2rc5  sin  (Pc5/2> 


r = 32. 7 in 
rC5 


q = 53. 6 deg 
C5 


i - 67  in 
*C5 


11  T 

<Or)  = (0,  -21.  S.  28) 


11P5 


cos  ri  r 
5 


0 sin  tc 

1 0 


P5.  PANEL  FIXED  FRAME 


-Sin  -1  r 
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-£ 
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Nominally, 


The  orienlalion,  location,  and  physical  dimensions  of  the  thruster  cone  surface  are 

defined  in  terms  of  the  parameters  (N_,,  v , v,  , O , h)  as  defined  in  Section  V.  In 

— e.  a b —a 

terms  of  B-frame  coordinates,  we  have 

( 

( 

+1 


> 


Y = 2.  6/2  = 15.0 
a 

= 3.2/2  = 18.5 
b 


B-Frame  origin  is  at  center  of  open  end 
of  thruster  cone. 


h = 23.  88 


Antenna  Rods 


Antenna  rods  are  simplified  as  cylindrical  surfaces  defined  in  terms  of  the  parameters 
(C,  L,  y,  l)  as  described  in  Section  II.  Antenna  rods  Aj  through  AJ2  are  arranged 
symmetrically  about  an  axis  L through  a point  as  shown  in  Figure  11-3.  Hence. 


_B 

(Lc.)  = (0,  0.  11 1 

_B 

(Cc>  = (8.  27,  0,  61.38) 


T 


For  Aj  through  A^.;,  i.e. , i = 1,12, 


Y.  - 1 . 25  in. 


1.  = 18.6  in. 

l 

-B  '] 

(L.  ) (0,  0,  1) 

— l 

B 


. T 


B 


(C.)  - (X.,  V.,  61.38)'  + (G  ) 
—l  ll  — G 


with 


X.  = -r  cos  9. 

l o i 

Y = - r sin  9. 

i o l 

r = 17.25  for  i = 1,  ....  8 
o 

= 6.  0 for  i = 7, 12 


= (i-1 ) x 45"  + 22.  5 ’ f or  i = 1 8 

= (i-9)  x 90"  + 45"  for  i = 9,  ....  12 


For 


r13  = 1-25in< 


= 44.3 
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Stowed  Vehicle  as  Cylindrical  Surface 


During  the  initial  spinning  phase,  the  GPS  configuration  has  its  solar  arrays  stowed. 

A cylindrical  surface  approximation  can  be  used  for  evaluating  solar  force.  The  cylin- 
der is  defined  in  terms  of  the  following  parameters: 

(L  ) = (0,  0,  1) 

_R  T 

(C  ) = (0,  0,  it 2) 

Y = 35.  7 

i = 67 


eg  Locations 

To  evaluate  the  torque  due  to  solar  pressure,  the  moment  arm  is  the  vector  differ- 
ence of  location  of  center  of  force  and  the  location  of  vehicle  center  of  mass.  In 
B-frame  coordinates,  the  locations  of  center  of  mass  for  GPS  are  summarized  in 
Table  11-1. 

O the r Si mplifi cations 

The  approximating  surfaces  are  assumed  to  have  uniform  physical  characteristics. 
For  GPS  spacecraft  exterior  surfaces,  the  coating  material  is  of  the  Al-Teflon  type 
which,  for  solar  arrays,  consists  of  a cover  of  AR  coat  and  a layer  of  UV  filter  coat. 


TABLE  H-l . SEQUENCE  MASS  PROPERTIES 


Principal  axis  angular  offset 


APPENDIX  I 


SOI. AH  ARRAY  FLEXIBILITY  MODEL 


l 
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APPENDIX  I 


SOLAR  ARRAY  FLEXIBILITY  MODEL 


A STARDYNE  (Control  Data  version  of  NASTRAN)  finite  element  program  was  used  to 
generate  a mathematical  model  of  the +Y  solar  array  model.  Only  the +Y  array  was  modeled 
with  the  pivot  point  considered  to  be  fixed.  The +Y  array  is  structurally  identical  to  the 
-Y  array.  However,  since  each  array  has  indeoendent  pitch  drive  actuators,  the 
model  considers  each  array  to  be  responding  separately  and  independently.  The  model 
will  be  discussed  first,  with  respect  to  the  original  coarse  model  that  was  used  for 
determining  the  model  structure  (i.c.,  how  many  structural  modes  should  be  included), 
followed  by  the  final  fine  mesh  model. 

COARSE  MESII  SOLAR  ARRAY  FLEXIBILITY  MODEL 

The  coarse  NDS  solar  array  structural  model  is  shown  in  Figure  1-1.  Nodes  5 through 
8 arc  on  the  boom  centerline,  at  the  approximate  location  of  the  solar  panel  to  boom 
hinges.  The  solar  panel  sun  sensors  are  located  close  to  Node  8. 


DEPLOYMENT  ARM  NODES 


BOOMS  : NODES 


©-KD 

©-<!) 


LUMPED  MASS,  INERTIA  OF  SOLAR  PANELS  NODES 


Firui  i T-i  Math  Model--Coarse  Mesh 


Figures  1-2  and  1-3  present  a layout  schematic  drawing  of  the +Y  array  and  the  Rockwell 
vibration  test  panel  orientation,  respectively.  Based  on  this  test,  the  equivalent  viscous 
damping  ratio  for  one  mode  was  estimated  to  be  3.  5 percent  (Zeta  = 0.  035).  This 
damping  estimate  was  obtained  informally  from  Rockwell. 


The  solution  for  solar  array  dynamic  response  can  be  obtained  by  direct  integration  of 
the  equations  of  motion  of  the  structural  math  model  shown  in  Figure  1-1.  However, 
such  a model  would  have  11  independent  nodes  with  six  degrees  of  freedom  per  node. 
Instead,  a solution  can  be  obtained  using  a reduced  set  of  generalized  coordinates.  The 
generalized  coordinates  are  the  natural  frequencies  (rad/sec)  of  the  structure.  The 
transformation  from  the  generalized  coordinates  to  the  physical  math  model  are  the 
normal  modes  (mode  shapes)  of  vibration. 

The  deployed  NDS  solar  array  flexibility  model  is  represented  by  a selected  set  of  n 
eigenvectors  and  eigenvalues.  The  solution  for  the  eigenvector/eigenvalue  problem  is 
based  on  a finite  clement  structural  model  of  the  NDS  deployed  solar  array  (+Y  array) 
based  on  available  structural  drawings  obtained  from  Rockwell  through  the  SAMSO 
contract  monitor  office. 

The  numerical  values  tabulated  in  this  memo  are  for  a coarse  mesh  math  model  based 
on  out-date  drawings.  The  associated  natural  frequencies  are  known  to  be  higher  than 
the  actual  solar  array  structure,  based  on  informal  results  of  a ground  vibration  test 
performed  at  Rockwell.  However  this  model  is  . ufficiently  close  in  dynamic  response 
characteristics  to  allow  checkout  r uns  to  be  made. 

EQUATIONS  OF  MOTION 

For  the  NDS  solar  array  model,  the  base  of  the  STARDYNE  elastic  structure  is  at  the 
following  station  (+Y  array): 

X = 0.  0 

V - 28.  77  in. 

Z = 38.  00  in. 

This  point  is  the  location  of  the  inner-  pivot  point  of  the  deployment  mechanism  and  is 
the  centerline  of  the  circular,  integrally  machined  drive  tube  as  shown  in  Rockwell 
drawing  No.  V505-932004. 


i 

i 


X 


TEST  NATURAL  FREQUENCIES 


(DEPLOYED +Y  ARRAY) 

f„,  Z AXIS  =4.1  Hz 

fn,  X AXIS  = 4.8  Hz 

f„,0Y  (TORSION)  =3.1  Hz 


Figure  1-3.  Rockwell  Test  Results’ 


For  convenience,  the  attach  point  of  th t'  elastic  structure  will  be  considered  to  be  on 
the  vehicle  centerline'  at  Z station  38.0.  Thus  modal  participation  factors  and  the  array 
response  torques,  moments,  and  forces  will  be  transformed  to  X = 0.0,  Y = 0.0,  and 
Z = 38. 00. 


This  will  simplify  the  equations  required  to  couple  the  elastic  array  model  to  the  rigid 
body  controls  simulation  model. 


i 


Reference:  Rockwell  report  N18231-26,  2/17/77. 


The  equations  of  motion  in  matrix  notation  for  the  +Y  array  are 


fm]  fii]  +fc]  fill  + fk]  [u]  = [m]fT]  [n] 


where 


are  motions  of  panel  points  on  the  array  relative  to  the  base  point.  The  accelerations 
of  the  base  point  are  [u].  { u]  has  six  possible  components:  X,  Y,  Z,  flX,  fi  Y,  az  . 

Let  the  panel  points  on  the  array  be  described  by  the  normal  mode  coordinate  trans- 
forma tion: 

{u]  = M {11]  (1-2) 

{ul  = M [il]  (1-3) 

fii]  = ftf]  fil}  (1-4) 

where  \&]  represents  the  normal  modes  and  the  {tl}  are  a set  of  generalized  coordinates 
which  are  uncoupled,  second-order,  ordinary  differential  equations. 

The  response  of  the  n^  "mode"  or  generalized  coordinate  is  expressed  as 


T)  + 2C  D Tl  + 0)“T1  = - fv  3 { ul 

n n n n n n 1 n 


where 


= equivalent  viscous  damping  factor  assuming  3.5  percent  (£  = 0.035) 


do  = nm  natural  frequency  (rad /sec) 


fy  1 

n g - model  participation  factor  for  base  acceleration  in  each  of  the  six  input 
axes:  X,  Y,  Z,  0X.  QY , 9Z. 


i 

t 


Equation  (1-5)  must  be  integrated  with  the  satellite  main  body  equations  of  motion.  If 
each  array  is  represented  by  the  first  four  elastic  modes,  then  four  sets  of  equations 
such  as  Equation  (1-5)  must  be  integrated.  If  solar  arrays  are  pulsed  separately,  then 
a total  of  eight  such  equations  will  be  required.  The +Y  array  is  assumed  identical  to 
the -Y  array.  However,  the  participation  factors  will  be  different  (of  opposite  sign)  for 
rotational  base  accelerations  in  9X  and  @Z. 

The  rotational  position  of  the  sun  sensor  at  the  outer  end  of  the  array  can  be  obtained 
by  using  Equation  (1-2). 

ARRAY  REACTION  FORCES  AND  MOMENTS 


Array  dynamic  responses  are  coupled  to  the  spacecraft  base  point  by  the  elastic  forces 
produced  at  the  attach  point  by  array  relative  responses:  [ul  = 


For  each  generalized  coordinate  T),  a set  of  three  forces  and  three  moments  are  pro- 
duced at  Node  1 in  the  STARDYNE  math  model.  These  have  been  transformed  to  the 
spacecraft  centerline  at  Z station  38.0.  Thus,  for  the  +Y  array. 


F 


M 


M 


M 


7.  * 


-VA2 
-P 
V A3 
-MA2 
-Torque- 


base,  L v 

nth 

mode 


M A3 


♦ 


STARDYNE 
Ream  1 
Notation, 
n1*1  mode 


+28.  77  VA3 

0 

^ +28.  77  VA2  J 


Modal  reactions  at  base  point 
per  unit  model  deflection  of 
+Y  elastic  array 
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For  the  -V  array,  the  modal  (or  generalized)  reaction  forces  at  the  base  point  (X  = 0.  0 
Y =0.0,  Z = 38.0)  are  as  follows: 


f > 

F 

X 

" 

r a 

-VA2 

0 

- 

F 

y 

+P 

0 

F 

VA3 

0 

Z 

= 

+ \ 

M 

X 

+ MA2 

-28.  77*VA3 

M 

y 

-Torque 

0 

M 

v ■/.  J 

-MA3 

v.  J 

+28.  77*VA2 

base  v ^ / 

Modal  reactions  at  base  point 
per  unit  modal  deflection  of 
-Y  elastic  array. 


SIGN  CON VKNTION 


The  sign  convention  for  the  loads  acting  on  spacecraft  and  the  loads  acting  at  +Y  array 
drive  point  is  shown  in  Figure  1-4. 


Figure  1-4.  Structural  Analysis  Sign  Convention 


6 


For  the  +Y  array  (STARDYNE,  Beam  1 loads,  end  A), 


F = -VA2, 

X 

+ XI  tip  deflection  gives  +F^ 

Fy  " -P' 

+ X2  deflection  gives  +F^ 

F = VA3, 
z 

+X3  tip  deflection  gives  +F^ 

M = -MA2, 

X 

+X3  tip  deflection  gives  +Mx 

M = -Torque, 

y 

+0Y  deflection  gives  +M^ 

M = M A3, 

7. 

+X3  tip  deflection  gives  -M 

For  the  -Y  array, 


F = VA2 

X 

+X1 

tip  deflection  gives  +F^_ 

Ph 

II 

^5 

+X2 

deflection  gives  +F 

y 

F = VA3 

7. 

+X3 

tip  deflection  gives  +F 

M = MA2 

X 

4X3 

tip  deflection  gives  +M^ 

M = -Torque 

y 

M = -M  A3 
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EQUATIONS  OF  MOTION,  FLEXIBLE  APPENDAGES 
COUPLED  TO  A RIGID  SPACECRAFT 


The  previous  derivation  may  now  be  extended  to  a general  math  model  as  follows. 


The  total  displacement  set  for  the  spacecraft  is 


where  the  spacecraft  rotations  [0^]  are  a subset  of  f U ^ 1 . For  no  external  loads  acting 
on  structure,  the  undamped  equations  of  motion  for  displacement  set  {U  ] become 

fM]fu‘)  + |'K]{u'}=0  (1-6) 

aa  a aa  a 
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Equation  1-6  may  he  partitioned  as  follows: 


(1-7) 


where  the  subset  {L'rl  of  the  displacement  vector  is  sufficient  to  define  rigid  body 
motion  without  additional  constraints. 


Now  introduce  the  following  "hybrid"  transformation  between  modal  coordinates  {11} 
and  physical  coordinates  [u]: 


Note  that 


K 

( relative 
to  base 


rVtV 


(I-8a) 


(I-8b) 


where 


[^  ] = matrix  of  eigenvectors  of  structure  with  { L^}  = 0 

|n  I = rigid  bodv  transformation  matrix 

i jjrJ 


= fT]  in  STAR  DYNE  write-up 


fl]  = identity  matrix 

^ } = ff  ] = vector  of  flexible  body  modal  (or  generalized)  coordinates. 


% 
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where 


7 


M(x)  (t),  My(t),  and  M_yt)  are  jet  torques  or  other  external  torque  applied  to 

spacecraft. 


NUMERICAL  VALUES  OF  FLEXIBILITY  COEFFICIENTS 


Table  1-1  presents  the  equations  of  motion  in  terms  of  generalized  coordinates  which 
must  be  integrated  for  each  solar  array.  The  coefficients  for  the  first  four  modes  are 
presented.  For  checkout  runs,  only  the  first  three  modes  are  necessary. 


TABLE  1-1.  NUMERICAL  VALUES  FOR  GENERALIZED  COORDINATE 
SOLUTION,  NDS  COARSE  MESH  SOLAR  ARRAYS* 

a.  +Y  Array 

[ I ]f  n}  + [ 2Cwn]{f|]  + [ w“  ] {T|}  = -f  I1  [u] 


f 1 ] = unity  matrix;  C = 0.03!!  or  less 


Mode 

No. 

— 

!•’ 

n 

(Hz) 

w 

n 

(rad /sec) 

+Y 

Array  Participation  Factors  (Vi)  for  Each  Mode 

X 

Y 

7. 

11  ♦ 

-X 

11  + 

6 V 

n+ 

ez 

1 

3.  074 

m.:m 

0 

0. 1676 

1.501 

10. 982 

0 

0 

2 

3.  491 

21.933 

1.  245 

0 

0 

0 

-0. 7963 

-9.  59 

3 

4.  334 

28.  491 

0. 0337 

0 

0 

0 

0. 4983 

-0. 344 

4 

10.  05 

62. 171 

0. 0483 

0 

0 

0 

-0. 2149 

-0.  155 

b.  -Y  Array 


-A 

.Array  Participation  Factors  (yi)  Tor  10a ch  Mode 

X 

A' 

7. 

ex 

e v 

a7 

HI 

i 

-10.  983 

0 

0 

2 

0 

-0. 7963 

+9.  59 

for  +\  Arr. 

3 

0 

0. 4983 

0.  344 

4 

0 

-0.  2149 

0.  155 

Reference:  Equation  (I- 1 0). 

STARDYNE  DYNRE2  values  divided  by  12.0  to  convert  inches  to  feet. 
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Table  1-2  presents  the  model  displacements  at  node  point  8 in  the  STARDYNE  model 
so  that  the  sun  sensor  displacements  and  angles  relative  to  the  base  point  can  be 
computed.  If  only  sun  sensor  rotations  are  needed,  use  only  the  values  in  the  columns 
titled  0X,  0Y,  and  0Z. 


TABLE  1-2.  NUMERICAL  VALUES  FOR  SUN  SENSOR  DISPLACEMENTS 

AND  ROTATIONS  PER  UNIT  GENERALIZED  DISPLACEMENT" 


a.  +Y  Array 


Pi 

Modal  Displacements  at  Node  8,  Coarse  Model 

e 

F 

n 

X 

(ft /ft) 

Y 

(ft/ft) 

Z 

(ft /ft) 

ex 

(rad /ft) 

0Y 

(rad/ft) 

ez 

(rad/ft) 

1 

3.  074 

0 

0. 0212 

1. 00 

0.  0209 

0 

0 

2 

3.  491 

0.  9825 

0 

0 

0 

-0. 0088 

-0. 0169 

3 

4.534 

1. 000 

0 

0 

0 

0.  2087 

-0.  0198 

4 

10.  05 

-0. 1109 

0 

0 

0 

0. 1284 

0.  0125 

Modal  Displacements  at  Node  8,  Coarse  Model 


9X 

(rad/ft) 

9Y 

(rad/ft) 

ez 

(rad  / ft) 

-0. 0209 

0 

0 

0 

-0. 0088 

+0. 0169 

0 

0.  2087 

+0.  0198 

0 

0.  1284 

0.  0125 

Reference:  Equation  11-7). 

f Example:  0Y  . =[0.  -0.0088  0.2087  0.1284] 

+v  array  1 

where 

T],  = generalized  displacement, 

1 i^1  mode 
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Tables  1-3  and  1-4  present  the  forces  and  moments  applied  to  a point  at  the  intersection 
of  the  solar  array  pitch  drive  axis  and  the  spacecraft  centerline  due  to  unit  generalized 
deflections  of  each  solar  array.  These  forces  and  moments  should  be  resolved  into  the 
torque  summation  about  the  spacecraft  eg.  For  use  with  elastic  arrays,  the  rigid-body 
mass  of  the  spacecraft  and  its  eg  should  not  include  the  solar  arrays. 

TABLE  I-3a.  COARSE  MESIT  MODEL  INTERNAL  FORCES  AND  MOMENTS,  BEAM  1, 
NODE  1 FOR  FIRST  FOUR  MODES  (STARDYNE)" 


Beam  Load 

Internal  Beam  Loads 

Conversion  to  Loads 

Component 

Mode  1 

Mode  2 

Mode  3 

Mode  4 

Acting  on  Base  Point 

VA2,  lb 

0 

-26. 619 

-19.044 

-155. 82 

FX  = -VA2 

P,  lb 

-1.  96 

0 

0 

0 

1 

FY  = -P 

VA3,  lb 

17.  5628 

0 

0 

0 

FZ  = VA3 

MA2,  in-lb 

-1036,  8 

0 

0 

0 

MX  = (-MA2  +28.  70  VA3) 

Torque, 

in-lb 

0 

+204. 335 

-3380. 29 

8319. 51 

MY  = -Torque  ^ 

MA3,  in-lb 

0 

-1695.  2 

-1787. 76 

-1530.01 

MZ  =^<MA3  +28.77  VA2) 

TABLE  I-3b.  FORCES  AND  MOMENTS  AT  (X  = 0,  Y =0,  Z = 38.  0 inches) 
DUE  TO  UNIT  GENERALIZED  DISPLACEMENTS 


Force 

Component 

Units 

Internal  Beam  Loads 

Mode  1 

Mode  2 

Mode  3 

Mode  4 

FX 

lb 

0 

26. 619 

19. 044 

155. 82 

FY 

lb 

1. 96 

0 

0 

0 

FZ 

lb 

1 7. 5628 

0 

0 

0 

MX 

ft -lb 

128.  507 

0 

0 

0 

MY 

ft -lb 

0 

-17.028 

281.69 

-693. 29 

MZ 

ft-lb 

0 

-205. 09 

-194. 64 

-501.08 

Reference:  Run  SSNSOUC 
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TABLE  1-4  FORCES  AND  MOMENTS  AT  X = 0,  Y = 0,  Z = 38.  0 INCHES 

DUE  TO  UNIT  GENERALIZED  DISPLACEMENTS  OF  -y  ARRAY 


Force 

Component 

Units 

Mode  1 

Mode  2 

Mode  3 

Mode  4 

FX 

lb 

0 

26. 619 

19. 044 

155. 82 

FY 

lb 

-1.96 

0 

0 

0 

FZ 

lb 

17. 5628 

0 

0 

0 

MX 

ft-lb 

-128. 507 

0 

0 

0 

MY 

ft-lb 

0 

-17. 028 

281.69 

-693. 29 

MZ 

ft-lb 

0 

205. 09 

194. 64 

501. 08 

If  the  solar  array  Z axis  is  shifted  more  than  5 or  10  deg  from  the  spacecraft  Z axis 
due  to  pitch  rotation  of  the  arrays,  then  the  appropriate  coordinate  transformations  are 
required.  In  fable  I - 1 , the  array  participation  factors  for  X,  Z,  9Z,  and  CY  must  be 
rotated.  Also  the  FX,  FZ,  MX,  and  MZ  forces  and  moments  in  Tables  I-3b  and  1-4 
must  be  rotated. 


CHECKOUT  RUNS 

Figures  1-5  through  1-7  represent  transient  responses  obtained  at  node  8 due  to  transient 
angular  acceleration  pulses. 

For  Figure  1-5,  the  array  pitch  torque  is  assumed  to  be  1.0  ft-lb  about  9Y  for  50  msec. 

< 2 

Using  the  Rockwell  I inertia  value  of  2.  2335  slug/ft  , this  causes  an  angular  accel- 

yy  2 2 

eration  of  0.4477  rad/sec“.  A reverse  input  of  -0. 153846  rad/sec  was  applied  for 


325  msec  to  kill  off  the  rigid  body  velocity, 
array  drive  pulse  input. 


Phis  forcing  function  approximates  a solar 


Two  STARDYNE  simulations  were  made,  one  using  the  first  three  modes  (up  to  4.5  Hz) 
and  a second  run  using  the  first  four  modes  (up  to  10.0  Hz). 
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DUE  TO  DYNAMIC 
AMPLIFICATION 


Figure  1-5.  Relative  Angular  Displacement  in  Torsion  of  NDS  Coarse  Model 
to  a Single  Drive  Torque  Pulse,  STAR  DYNE  DYNRE1  Solution 


X x 102.  in 


4486  ft  lb  AT  t - 120  uc 


Figure  1-6.  Relative  Room  Tip  X Displacement,  NDS  Coarse  Model, 
Due  to  Thruster  Jet  Pulse,  STARDYNE  DYNRE1  Solution 


Reference:  Run  SSNSOQ9 


123 


0y  x 104,  rad 


(P° 

O O 


RESPONSE  FOR  FIRST  THREE  MODES 
f = .035 


.0053  rad/seir 


.3  o -4 


>-  .10 


CL  O 

cbc P° 


MAX  ROTATION  

= -.2132  x 10'4  rad 

AT  t = .42  sec 


MAX  TORQUE  REACTION  TO 
BASE  (+  Y ARRAY) 

= -.449  in-lb  = -.0374  »t-lb 
AT  t = 42  sec 


Figure  1-7.  Room  Tip  and  Sun  Sensor  9Y  Rotation  Due  to  Yaw  Thruster 
Jot  Pulse.  STAR  DYNE  DYNRE1  Solution" 


Figure  1-6  presents  the  simulation  for  the  arrays  excited  by  the  thruster  jet  firing  about 

the  yaw  axis.  Using  an  inertia  estimate  of  100  slug  ft"  and  a torque  of  0.53  ft-lb,  an 

2 

acceleration  ot  0.005.1  rad/sec  about  the  array  Z axis  was  input  into  the  array  base. 
The  X response  of  node  8 at  the  tip  of  the  boom  is  plotted. 


Figure  1-7,  also  for  the  same  yaw  thruster  pulse,  shows  the  +y  array  torsional  response 
which  is  coupled  with  the  X deflections  due  to  the  deployed  boom  offset  from  the  pivot 
axis.  The  DY  NR  El  solution  was  carried  out  to  2.5  sec  and  the  peak  ?Y  response  plotted 
at  0.  42  sec  is  the  maximum  displacement  observed. 


For  a roll  thruster  pulse  of  0.0834  ft-lb  with  a rigid  body  inertia  of  90  slug  ft  , Figure 
1-8  shows  the  transient  boom  tip  displacement  in  the  Z direction. 


Reference:  Run  SSNSOQ9 


124 


r 


Z„103.m  Z DIRECTION 


MAXIMUM  + V ARRAY  REACTIONS  ON  SPACECRAFT 
AT  X=Q.  YrO..  Z=38.0 

FZ  = VA3  = -.00899  lb  AT  t = .13  sec 
MX  = 1/12  (28.77  x VA3  - MA2) 

.0658  ft  lb  at  t = . 13  sec 

Figure  I-tt.  Relative  Boom  Tip  Displacement,  Z,  Due  to 
Single  Roll  Thruster  Pulse 


' Note  that,  these  checkout  runs  were  made  with  the  array  Z axis  oriented  along  the  space- 

I craft  Z axis. 

FINE  MESH  SOLAR  ARRAY  FLEXIBILITY  MODE 

After  the  solar  array  flexibility  structure  had  been  determined  and  the  latest  boom  and 
array  stiff  data  obtained,  a fine  mesh  solar  array  flexibility  model  was  constructed. 
Table  1-5  presents  numerical  values  for  the  generalized  or  modal  coordinates  of  the  +Y 
and  the  -\  array.  (Separate  array  models  are  used  since  they  have  individual  pitch 
drive  control  systems  and  are  not  physically  connected  for  torsion.) 

Table  1-6  presents  the  appropriate  mode  shape  displacements  and  rotations  at  the  outer 
ends  of  each  boom  (Y  = 104.97  inches)  which  is  close  to  the  sun  sensor. 


Table  1-7  presents  the  coupling  terms  which  are  used  in  torque  contributions  due  to 
array  generalized  accelerations. 


n 

TABLE  1-5.  NUMERICAL  VALUES,  GENERALIZED  COORDINATE  SOLUTION, 
NDS  FINE  MESH  SOLAR  ARRAYS* 


a.  +Y  Array 

[ini’]  +[2Cwn][Tl]  + r ]{T1}  = -[  ]T[n) 

[ I ] = unity  matrix;  Q = 0.  035  or  less 


Mode 

No. 

F 

n 

(Hz) 

W 

n 

(rad/sec) 

X 

1 

3.9133 

24. 588 

0 

2 

3. 9138 

24. 591 

1.3305 

3 

5.  086 

31. 958 

0. 3143 

4 

9.  75 

61. 289 

0 

Participation  Factors,  fvl  , for  Each  Mode 


Z 


ex* 


0Y 


ezf 


0.  2140 
0 
0 

-0. 2788 


1. 8102 
0 
0 

0.  4589 


12.  159 
0. 00039 
0 

0. 1847 


0.  00005 
-0.  95706 
2.  6012 
0 


0.  00048 
-9. 5772 
-2. 8584 
0 


b.  -Y  Array 


Mode 

F 

W 

Participation 

T 

Factors,  fy)  , for  Each  Mode 

No. 

(Hz) 

(rad/sec) 

X 

Y 

.. 

Z 

ex 

8Y 

ez 

1 

-12.  159 

0. 00005 

2 

or  +y  Arraj 

-0.00039 

-0. 95706 

9.  5772 

3 

0 

2.6012 

2. 8584 

4 

-0. 1847 

0 

0 

* 

Reference:  Equation  (1-5). 

*STARDYNE  DYNREI2  values  divided  by  12.0  to  convert  inches  to  feet. 
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TABLE  1-6.  FINE  MESH  SOLAR  ARRAY  MODEL.  MODE  SHAPE  AT  SUN  SENSORS 
(Displacement  and  Rotations  per  Unit  Generalized  Displacement) 


a.  +Y  Array 


Mode 

F 

Modal  Values  at  Node  8,  Fine  Mesh  Model 

Number 

n 

X 

(ft /ft) 

Y 

(ft /ft) 

Z 

(ft /ft) 

0X 

(rad/ft) 

8Y 

(rad/ft) 

ez 

(rad /ft) 

1 

3. 9133 

-0. 00031 

0. 03157 

0. 7797 

0.01224 

0 

0 

2 

3. 9138 

0. 9613 

0 

0 

0 

-0. 00697 

-0.01453 

3 

5.  086 

0. 2432 

0 

0 

0 

0. 01852 

-0. 00388 

4 

9.  75 

0 

-0. 0172 

-----  - 

-0. 437 

-0. 00709 

0 

0 

b. 


-Y  Array^ 


Mode 

Number 

F 

n 

Modal  Values  at  Node  8,  Fine  Mesh  Model 

X 

(ft /ft) 

Y 

(ft /ft) 

Z 

(ft /ft) 

ex 

(rad /ft) 

0Y 

(rad /ft) 

9Z 

(rad/ft) 

1 

-0. 01224 

0 

0 

2 

0 

-0. 00697 

0. 01453 

— 

Same  as  +Y  An 

•ay  (above)- 



3 

0 

0. 01852 

0. 00388 

4 

0.  00709 

0 

0 

.Reference:  Equation  (1-2). 
’Reference:  Run  SSNSOQY 
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TABLE  I-7a.  NDS  FINE  MESH  SOLAR  ARRAY  MODEL  COUPLING  TRANSFORMATIONS 
FOR  SPACECRAFT  TORQUE  SUMMATIONS  DUE  TO  GENERALIZED 
ACCELERATIONS,  FIRST  FOUR  MODES 


Mode 

No. 

M 

n 

General 

Mass 

(lb-sec^/in) 

f Mn  ]|Y 
n m 

] - f M ] (Participation  Factors) 

X 

Y ' 

z ■ 

9X 

BY 

BZ 

1 

0. 028288 

0 

0.  00608 

0. 0512 

4. 1274 

0. 000016 

0. 00016 

0.044165 

0. 05876 

0 

0 

0.000208 

-0. 5072 

-5. 0756 

0. 022699 

0. 0071 

0 

0 

0 

0.7085 

-0. 7786 

0. 011422 

0 

1 

O 

o 

o 

[C 

0. 0052 

0. 0253 

0 

0 

— 

I Divide  by  12  if  used. 

TABLE  I -7b.  f Mn  ]f  y ] ' FOR  TORQUE  SUMS  DUE  TO 
i ' . GENERALIZED  ACCELERATIONS 


Axis 

Mode  1 

Mode  2 

— 

Mode  3 

Mode  4 

Array 

TXX 

4.  1274 

0.  000208 

0 

0. 0253 

TYY 

0.  000016 

-0. 5072 

0. 7085 

0 

TZZ 

0. 00016 

-5. 0756 

-0. 7786 

0 

TXX 

-4.  1274 

-0. 000208 

0 

-0. 0253 

TYY 

+0. 000016 

-0. 5072 

0. 7085 

0 

-Y 

TZZ 

-0.  00016 

5 . 07  5 6 

0. 7786 

0 
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Figures  1-9  through  I- 1 1 present  geometry  plots  of  the  fine  mesh  math  model.  The 
deployment  arm  and  fittings  were  not  changed  from  the  previous  model  (see  Figure  1-12). 
The  honeycomb  aluminum  solar  panels  were  modeled  as  shown  in  Figure  1-13. 

Triangular  plate  elements  were  used  because  the  STARDYNE  quadrilateral  plate  element 
cannot  model  a honeycomb  panel. 

Each  panel  has  a set  of  four  hinges,  similar  to  door  hinges,  to  allow  folding  ior  launch. 
At  each  end  of  the  panel  a linkage  is  used  to  hold  the  panels  out  in  the  deployed  position 
for  orbit.  These  links  were  modeled  by  scaling  off  the  drawings,  and  some  error  in 
the  torsion  mode  could  be  caused  by  inaccuracies  in  the  linkage  model. 


6 


Figure  1-9.  NOS  +Y  Array  Fine  Mesh  Dynamic  Model,  X3  Rotated  View 
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XI  VIEW  ROTATED  AXES 

.35  22.10  31  .36  42.51  53.36  £+.12  14.86  85.61  96.35 


/-<  iK-ChJd  Uk  ( 4 pi„£M) 

/ CCO-^oifi  X U*,  > 

vjUi^  dtpU^c^ 


13.114  -Y-l.kl  - T1 5 -00  \,.bi,  ~T.1'1  ZZ  .03  I’z.OC 


Figure  I— 11.  YDS  +Y  Array  Mode]  Knd  View,  I,oo king  Down  Y Axis 


SUPPORTING  DATA 


Origin  for  Array  Cylindrical  Coordination  System 

ICYL  = 2 
X2  = H 

A = 0 deg  at  +X3  axis 

X3  = Z station  = 33.  0 

eg 

approximately  1.  25  in.  above  inside  mean  radius 

R *=32.25 

mean 


X3  = (32.  25  - 1.  25)  + 33.  0 


= 64. 

. 0 in. 

(Adjust 

origin 

after 

geometry 

check) 

For  -\ 

Array 

Point  A 

: XI  = 

-17.  7 

Measu 

red  Drawing  No. 

. V505- 

932003 

X2  = 

0 

From 

angle 

=180  deg 

+ 28.  32 

= 208. 

32 

X3  = 

64. 

To  angle  = 1 

80  deg  - 2 

8.32  = 

151.  68 

Point  if 

: XI  = 

-0.  0 

Point 

C:  XI 

= -10. 

X2  = 

0.  0 

X2 

= +10. 

X3  = 

64. 

X3 

= 64. 

For  +X 

Array 

Point  A 

: XI  = 

+ 17.  7 

Point 

B:  XI 

= 50. 

Point  C 

::  xi  = 

30. 

X2  = 

0. 

X2 

= 0. 

X2  = 

10. 

X3  = 

64. 

X3 

= 64. 

X3  = 

64. 

TABLE  1-8.  ROCKWELL  WEIGHTS  DATA 


x-r. 


One  Solar  Panel  (-X  Array) 
Weight  = 15.  34  lb  excluding  boom 


X = -15.  5 = XI 

eg 

Y = 71.3=  X2 

eg 

Z = 33.0  =X3 

eg 


) STAR  DYNE  Notation 


Total  area,  one  panel  = 2133.  576  in~ 

Modify  material  density  (MATLG  2 card)  to  make  face  sheet  density  match  total  density, 

Total  volume,  face  sheets  = AREA  x 0.020  in. 

= 2133. 576  x 0. 020 
= 42. 67152  in3 


_ . 15.34  1b 

Composite  gross  density  = 

p9  = 0.  35949  lb /in3 


Solar  Panel  Distributed  Weights 


Honeycomb  core  material  = 2.3  lb / ft 


D : T~  0.001331  lb/in" 


Composite  density  based  on  face  sheet  thickness  is,  for  1 square  inch  of  panel. 


wt-  =0.10  lb/in'  0.  020  in  = 0.  002  lb 
face 

wt  = 0.001331  x 0.50  = •^-°^;)plb 
core  .0026655 
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Increase  density  of  face  sheets; 


Q 


0. 0026655 
0.002 


x 0.  10  = 0. 13328  lb/in 


3 


to  account  for  core  material  and  face  sheets. 


Distributed  Weight,  One  Panel 

2 2 

AREA  = 72  tri-plates  x 29.  633  in  = 2133.  576  in~ 


From  weights  statement,  -X  panel 


GPS  FAA  ARRAY  = 6.92  lb 

FA  A 450002  SUBSTRATE  = 8.  22  lb 

FDB  SLR  PNL  SUBSTRATE  = y--||  lb 

Weight  Shifts  on  Solar  Panel  (STAR DYNE  Run  SSNS05A) 


(-X  panel  geometry  check) 


WT  = 15.  34  lb 


X 


eg 


Y 


eg 


v Z 
v eg 


-17. 7 in. 
71.47  in. 
33. 1173  in. 


based  on  uniform 
density  panels 


Desired  eg: 


X 

eg 

Y 

eg 


Z 


eg 


-15.  5 in. 
71.3  in. 
33.  0 in. 


The  model  can  be  simulated,  from  the  distributed  weights  above,  as  follows 
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Node  No. 


Weight 

(lb) 


X,  in. 


21 

41 

61 

81 

101 

121 

141 


0. 10653 
0. 21306 


-33.0 


▼ 

0. 10653 


W • X 

-3. 51549 
-7.03098 


▼ 

-3. 51549 


Totals : 


1. 27836 

= W . 
edge 


W • X = 15.  34  x -17.  7 


-42. 18588 


-271. 518 


- current 


desired  W • X = 15.  34  x -15.  5 = -237.  77 


MW  ■ X)  = 33.  748 


W , x (XI  - XI)  = 33.748  in-lb 
edge  new 

(XI  - (-33.  0))  = :~^7r  = 26.  3994  in. 
new  1 . 2783b 


XI  = 26. 39945  - 33. 0 = -6. 6005  in. 
new 


Weight  Shift,  -X  Panel  (X-Direction  Shift) 


W . = 1.27836  lb 

edge 


X , =-33.0  in. 

edge 


W . X - 15.34  x -17.7  = -271.518  (current) 

panel  c 

W *X  , = 15.34  x -15.5  = -237.77  (desired) 

panel  desired  


AW  • X = 33.748  (to  match  Rockwell  weight  report) 


W = a • W , = weight  left  at  outer  edge  of  panel 

1 edge 


W 


„ = (1  - a)W  , = weight  shifted  to  inner  edge  of  panel 

2 edge 


W (X 


?dge 
X = AW  • X 


new  old 
W.,(-2.  400744  + 32.  99926)  = 33.  748 

W.,(30.  59852)  = 33.  748  in-lb 


W = 1. 102929  lb 
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Subtract  1.  102929  lb  from  outer  edge. 
Add  1.  102929  lb  to  inner  edge. 


Node 

Node 

AW 

No. 

No. 

27 

31 

47 

51 

■ 

07 

71 

87 

91 

107 

111 

127 

131 

1 

r 

147 

151 

+0.09191075 

KI 

IT 

1 

| 

Node 

No. 

Node 

No. 

A 

W 

21 

37 

-0. 09191075 

41 

57 

-0. 1838215 
v y 

61 

77 

81 

97 

101 

117 

121 

137 

i 

141 

157 

-0. 09191075 

-1. 102929 

Aluminum  Tube,  Strut  Assembly  (Deployment  Arm) 


Major  portion  is  machined  down  to 


I ■ f H t 
x-x 


= tt  (2.  54/2)'*  0.040 
4 

= 0. 2574  in 


I I 0.2574 

7.-7.  X-X 

9 

A = rDt  = 0.  21 4 in“ 

I = ZI  0.  52  in4 
yy  xx 


BPROP1  No.  1 


ID  2.  500  in.  diameter 


t ,,  - 0.  040  in. 
wall 

mean  diameter  = ID  + t 


wall 
- 2.  540  in. 


Part  No.  V505-.340052,  Reference:  D/N  V505-932004. 
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k 

i 


■> 


I = 0.014  - ignores  2 x K • d“  stiffness 
A = (1.  50)(0.  25)(2)  = 0.  75  in" 


9 

Area  = 0.  0672  in'" 

X2  = ~r!|p-  = 0.  24125  1X2  = ^ 0.  63  x 0.  503  - ~ 0.  59(0.  42)3 

1X3  = 0. 0010965  + 0. 0016692  = 0. 0065625  - 0. 00364266 

=0.002766  in4  1X2  = 0. 00292  in4 

J = r BLt3  « i{0.  63  + 0.  63  + 0.  42)(0.  040)3 
= 0. 0000358  in4 

BPROP1  No.  11,  Reference:  Rockwell  Drawing  No.  V505-450003. 
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EFFECT  OF  DRIVE  LINK  ADDITION 


Because  of  the  discrepancy  in  model  versus  test  torsional  frequency,  the  effect  of 
adding  a flexible  torsion  lock  or  drive  link  was  investigated  using  the  coarse  mesh 
model.  I-'igure  1-14  shows  the  modified  model  with  the  flexible  drive  link  introduced 
[ with  nodes  28  and  29.  1 

Table  1-9  presents  the  results  of  such  modification  on  natural  frequencies.  The 
torsional  frequency  was  reduced  from  6.  84  Hz  to  from  5.  13  to  6.  16  Hz,  depending  on 
the  model  of  the  drive  link  used.  Although  test  and  model  frequencies  were  still  in 
slight  disagreement,  these  results  indicate  that  the  torsional  frequency  would  be  higher 
in  flight  than  the  test  indicates. 


I’ ABLE  1-9.  STAR DY NE  FREQUENCIES  OF  MODIFIED  MODFI. 


Mode  Shape 

Test 

STAR  DYNE  Mode 

Frequencies  (ID) 

Description 

Frequency 

(ID.) 

\ 

■ 

B 

C 

D 

Bonding, 

normal  to  panel 
Z-direction 

4.  1 

3.  07 

8.  75 

3.75 

3.75 

Bending,  in 

plane  to  panels 
\ -di  roc  t ion 

4.  8 

4 S' 

3.  75 

3.  70 

3.  74 

Torsion  about 
-axis 

8.  1 

4.  53 

6 . 84 

5 . 1 3 

r>.  1G 

C ondition: 

\ From  reference  (old  boom  stiffnesses). 
E With  new  boom  stiffnesses. 

C With  drive-link  addition. 

D With  expanded  drive-link  addition. 
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COMBINED  EARTH  SENSOR  MODEL 


APPENDIX  J 


COMBINED  EARTH  SENSOR  MODEL 


INTRODUCTION 

This  appendix  describes  the  model  developed  for  the  GPS  combined  earth  sensor. 
The  model  can  be  used  to  evaluate  acquisition  as  well  as  normal  operation. 

This  model  is  based  on  information  obtained  from  Barnes  Engineering  Company 
at  a meeting  in  St.  Petersburg,  Florida  on  March  13,  1975,  and  a subsequent 
phone  call,  and  the  Rockwell  procurement  specification  for  the  sensor. 

SENSOR  DESCRIPTION 

The  Combined  Earth  Sensor  is  a two- axis  , static  sensor  which  incorporates 
two  horizon  crossing  indicators  for  spinning  mode  attitude  sensing.  The 
adjective  static  means  the  sensor  uses  no  moving  parts  in  the  measurement 
of  attitude. 

The  two-axis  sensor  consists  of  an  optical  system,  24  thermopile  detectors, 
and  processing,  electronics.  The  24  detectors  are  arranged  on  the  focal 
surface  of  the  optical  system.  The  earth  image  is  tangent  to  the  inner 
edges  of  the  A detectors  at  null  at  the  design  altitude.  Two  roll  and 
two  pitch  measurements  are  computed  sequentially  by  switching  sums  and 
differences  of  various  detector  outputs  to  the  input  of  an  integrator. 

Each  measurement  is  sampled  and  held  for  comparison  check  and  possible 
output.  A new  roll  and  new  pitch  measurement  is  computed  each  1/4  second. 

Roll  and  pitch  computations  are  separated  by  1/S  second.  The  S detectors 
are  for  space  reference  and  arc  used  with  the  A detectors  for  accurate 
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measurements  over  +1°  in  the  track  mode.  The  A and  B detectors  are  used 
in  the  acquisition  mode  for  +4°  range.  The  thermopile  detectors  have  a 
time  constant  of  1.4  seconds  and  are  sensitive  to  the  14  to  16  micron 
wavelengths  (C0?  band) . The  sensor  continuously  tests  all  signals  and 
rejects  spurious  or  incorrect  signals  due  to  sun,  moon,  or  other  anomalies. 

The  two  horizon  crossing  indicators  consist  of  3 degree  by  3 degree  square 
fields  of  view.  One  diagonal  of  each  square  is  parallel  to  the  scan  path. 
One  field  is  10°  from  the  -X  axis  toward  +Z  and  the  other  is  10°  from  -X 
toward  -Z  (Spin  rotation  is  about  the  Z axis).  The  detectors  are  pyro- 
electric with  insignificant  time  constant. 

MODEL  DESCRIPTION 


The  model  consists  of  an  initialization  and  two  normal  operation  sections. 

The  initialization  section  performs  one-time  computations  and  clears  outputs 
and  counters.  One  normal  operation  section  is  for  the  horizon  crossing 
indicators  used  in  the  spin  phase.  The  other  normal  operation  section  is 
for  two-axis  earth  sensing. 

Characteristics  of  the  horizon  crossing  indicator  model  are  as  follows: 

• Inputs  to  the  model  are  spacecraft  position  vector  from  center 

of  earth,  and  spacecraft  attitude  direction  cosine  matrix.  Both 
inputs  must  be  relative  to  the  same  reference  frame. 

• When  a sensor  line  of  sight  crosses  the  edge  of  the  earth,  the 
time  from  the  crossing  to  the  sample  time  is  output.  If  no 
crossing  occurs  in  the  sample  interval,  the  output  is  set  to 
-At . 

• If  a sensor  line  of  sight  crosses  the  earth  very  close  to  a tangent, 
there  may  not  be  an  output.  The  sample.  At,  should  be  small 


enough  so  that  this  loss  of  output  does  not  happen  when  attitude 
determination  is  critical.  Select  At  as  follows: 


e 


where 

W = spin  rate  in  rad/sec 

s 

= — W 
60  RPM 

W„t,„  = spin  rate  in  rev /min 

RPM 

h = satellite  altitude  in  n.  mi.(  100  _<  h _<  10897) 

r = earth  radius  in  n.  mi.  (3440) 

e 

i = distance  from  earth  center  to  chord  of  detector 

Path  in  n.  mi. 


e.g. 

let 


W 


RPM 


h 


r 

e 

l 


= 100 
= 10897 
= 3440 

= 3405.6  = 0.99  r 

e 


then  At  < 0.00326  seconds. 


o The  equations  used  in  the  flow  diagram  are  derived  below 


R X LOSi 

1*1 


SIN  (ANG2) 


R • LOS 

~WT~ 


-COS  (ANG2) 


ANG2  = TAN  1 

SIN  ( ANG2 ) 

COS  (ANG2) 

ANG2  = TAN_1 

| R X LOS | 

-IT  • LOS 
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ANGl  = SIN 


n 


-1 


r 

e 

TT 


= TAN-1 


re/|R| 


2-ih 


[1  - ( re / I R i ) ] 


= TAN 


-1 


[IR]  2 - re2]h 


ANG3  = ANG2  - ANGl 


Note  that  3440  + 100  ‘R|  10897  + 3440  n.  mi, 


13.8°  < ANGl  < 76.4°. 


Characteristics  of  the  two-axis  earth  sensor  model  are  as  follows: 

• The  sensor  output  is  updated  every  1/8  second. 

• The  sensor  model  has  four  states  which  are: 


Compute  Roll  1,  Output  Pitch  2 if  valid 
Compute  Pitch  1,  Output  Roll  1 if  valid 
Compute  Roll  2,  Output  Pitch  1 if  valid 
Compute  Pitch  2,  Output  Roll  2 if  valid 

\ 

The  state  is  advanced  every  1/8  second. 


State  1 
State  2 
State  3 
State  4 


Each  detector  area  is  divided  into  1°  x 1°  squares. 


The  area 


coverage  of  each  1°  x 1°  square  by  the  earth  is  computed  from  the 


distance  between  center  of  the  earth  circle  on  the  focal  surface 
and  the  center  of  the  square.  The  orientation  of  the  square  is 
not  considered. 


The  position  of  the  earth  image  on  the  focal  surface  (B&C)  is 
computed  from  true  attitude.  The  diameter  of  the  earth  image 
on  the  focal  surface  (RE)  is  computed  from  satellite  altitude. 
Measurements  on  the  focal  surface  are  in  equivalent  radians. 
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Each  detector  output  is  the  sum  of  the  coverages  of  1°  x 1° 


squai-es  within  the  detector. 


The  sensor  output  is  computed  from  the  detector  outputs  according 
to  the  functions  in  the  electronics. 


Functions  of  the  system  programmer  to  reject  spurious  or  incorrect 
signals  is  included. 

Subroutine  must  be  initialized  by  calling  with  MODE  = 0 prior 
to  operating  in  two- axis  mode. 

The  time  increment,  At,  should  be  selected  from  one  of  the 
following  to  synchronize  with  the  sensor  update  frequency: 

Recommended  At 


0.125 

0.0625 

0.025 

0.0125 

0.01 

<0.01 


Synchronization  error 
0.005  sec.  is 
insignificant 


EARTH  RADIANCE  VARIATION 


The  variation  of  radiance  over  earth  surface  for  the  14-16  p C0o  band  has 
been  calculated  in  Figure  J-l.  Plots  of  the  radiance  as  function  of 
latitude  and  time  of  year  for  various  altitudes  is  shown  in  Figure  j_2.  The 

model  for  CPS  simulation  is  deduced  from  these  results.  For  maximum  radiance 
variation  a January  situation  is  considered.  The  radiance  is  assumed  to  be 
a function  of  latitude  but  not  longitude.  The  radiance  variation  for  south- 
ern hemisphere  is  assumed  to  be  that  for  northern  hemisphere  six  months 
later.  The  radiance  variation  for  earth  disk  is  approximated  by  piecewise 
linear  segments  using  the  values  corresponding  to  20  Km  altitude  shown  in 
Figure  J-2.  The  radiance  values  and  the  normalized  variations  are  summarized 


J 
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in  Table  J-l.  The  normalized  variation  is  used  to  be  compatible  with  the 
existing  earth  sensor  simulation  that  does  not  include  the  details  of 
sensor  optics  and  filter  response.  The  radiance  for  the  15°  latitude  band 
is  taken  as  that  corresponding  to  a 220 °K  earth  which,  in  conjunction  with 
the  sensor  optics  and  filter  response,  provides  3.15  uw  when  the  "A"  detector 
FOV  is  fully  illuminated. 


TABLE  J-l.  EARTH  RADIANCE  VARIATION  WITH 
LATITUDE  AND  TIME  OF  YEAR 


LAT. 

MON. 

W/M2 

BADIAHCE  STER 

NORMALIZED 

VARIATION 

-75° 

JAN 

4. 13 

0.72 

-60° 

JAN 

4.63 

0.80 

1 

•P- 

Ln 

O 

JAN 

4.88 

0.85 

-30° 

JAN 

5.5 

0.96 

+15° 

JAN/JUL 

5.75 

1.00 

30° 

JUL 

5.94 

1.03 

45° 

JUL 

6.19 

1.08 

60° 

JUL 

6.44 

1.12 

75° 

JUL 

6.94 

1.21 
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Implementation  of  the  earth  radiance  variation  into  the  sensor  model  requires 
the  following  computations  which  are  based  on  the  geometry  given  in  Figure 


(a)  LOS-vector  for  detector  cell  element  centered  at  (X  , Y ) : 

a a 

“ -VXa2  + "a' 

X sin  U 

X = 

s U 

-Y  sin  U 
„ a 


Z = cos  U 
s 


Denote 


(L  S) 
—a 


(b)  Look  point  on  earth  surface  along  L : 

3 


O - tib  Lb1 


a,1)  ■ (Ev1)  + s O 


whe  re 


(R  ) = look  point  on  earth  surface  along  L . 

3 3 

(Ry1)  = vehicle  position  from  center  of  earth 
S = slant  range  from  Ry  to 

= "(lyVO  - '>e2-^vI)T(^VI)  + f(lvI)T(^aI):,2}iS 

Y = 3590 . 95  n . mi . 
e 


(d)  Scaling  factor  for  radiance  variation: 

V = F(B) 

where  F(6  ) is  the  piecewise  linear  function  defining  the 
radiance  variation  with  latitude;  F(3)  is  plotted  in  Figure  .7-3. 

(e)  The  value  of  AREA  is  modified  from  the  original  value  by  the 
scaling  factor,  i.e.: 

AREA-*- AREA  * V 

Figure  J-/|  presents  a comparison  of  the  GPS  simulation  radiance  model  and  a 
model  described  in  a Barnes  Engineering  document,  "Technical  Description 
Model  13-16901  and  Model  13-16903  Synchronous  Altitude  Static-Sensors", 

22  August,  1975.  Models  agree  reasonably  well. 
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SENSOR  ERROR  MODELING 


The  sensor  error  modeling  considered  here  includes  detector  and  pre-amp 
noises,  electronics  parameter  shift  error,  error  due  to  temperature 
gradient  across  detector  and  sensor  alignment  error.  Horizon  uncertainty 
errors,  though  not  a sensor  originated  error,  are  also  included  in  the  con- 
sideration. These  are  errors  that  exist  in  all  operational  situations. 

Errors  due  to  sun  or  moon  in  sensor  detector  field  of  view  and  errors 
resulted  from  variation  of  earth  disk  size  due  to  vehicle  orbit  eccentricity 
exist  only  for  a specific  operational  condition  and  are  not  included  here. 

The  detector  and  preamplifier  noise  has  been  calculated  to  be  0.379  uV  (3c) 
for  sensor  in  track  mode  (A  detectors  only)  and  0.464  pV  (3')  for  sensor 
in  acquisition  mode  (A+B  detectors)  with  a detector  scale  factor  of  3.15  pW / 
deg.  and  responsivity  of  20  V /W . These  noises  are  equivalent  to  0.0060  deg. 
(3c)  and  0.0.074  deg.  (3o)  respectively.  Use  the  output  scale  factors  of 
5 V/deg  and  1.25  V/deg  for  track  and  acquisition  modes  respectively.  The 
noise  at  output  then  becomes  0.03  V(3 c)  and  0.0093  V(3o).  The  noise  error 
will  be  simulated  as  white  noise. 

Electronics  parameter  shift  has  been  analyzed  by  Barnes,  Reference  J-2. 

The  values  at  initial  setting  and  at  end  of  life  have  been  calculated  to  be 
0.0034  deg.  (3a)  and  0.0219  deg.  (3a)  respectively.  These  correspond  to 
output  errors  of  0.017  V (3a)  and  0.1095  V (3c)  for  track  mode  and 
0.0043  V (3o)  and  0.0274  V (3o)  for  acquisition  mode.  The  error  due  to 
electronics  parameter  shift  will  be  simulated  as  bias  errors. 

The  temperature  gradient  considered  by  liames  amounts  to  a tenperature 
difference  of  0.00048°F  between  A and  S reference  junctions.  The  correspond- 
ing attitude  error  evaluated  for  temperature  of  70°C  is  0.0103°  or  0.0515  V 
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for  track  mode.  For  acquisition  mode,  the  temperature  change  of  0.0028  deg. 
between  detectors  is  calculated  to  cause  attitude  error  of  0.0029  deg.  which 
amounts  to  output  error  of  0.0035  V.  All  values  quoted  above  are  considered 
to  be  for  worst  case  and  will  be  treated  as  3a  values.  The  error  due  to 
temperature  gradient  will  be  simulated  as  bias  errors. 

The  achievable  static  sensor  alignment  accuracy  has  been  estimated  by  Barnes 
to  be  in  the  order  of  0.015  deg.  (3c).  The  equivalent  output  error  is 
0.075  V (3c)  for  track  mode  and  0.019  V (3c)  for  acquisition  mode.  The 
effect  of  alignment  error  will  be  simulated  as  bias  errors.  The  horizon 
uncertainty  error  has  been  treated  as  an  exponentially  correlated  noise 
with  standard  deviation  0.88  Km,  correlation  time  10  days  and  correlation 
distance  2500  N.  Mi  in  Reference J-l,  i.e., 

E[h(t,s)h(thT,  S+d)]  = chVT/1°  e_d/250° 

where 

o,  = 0.88  Km 
T = time  in  days 
d = distance  in  n.  mi. 

-T/10 

The  term  e can  be  approximated  by  1 if  T is  much  less  than  10  days. 

Factoring  in  the  GPS  orbit  alticude  of  10897  N.  Mi  the  horizon  uncertainty 
error  amounts  to  0.0020  deg.  of  earth  disk  angle  error.  In  the  evaluation 
of  and  R^  detectors  at  opposite  sides  of  the  earth  disk  are  involved. 

The  separation  between  the  look  points  of  these  detectors  is: 

x 3461 .6  N . Mi  . 


(180°  - 14°) X- 

d'~ 7? 


= 10029  N.  Mi . 
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The  separation  is  approximately  four  times  the  correlation  distance. 

Hence  the  effects  of  horizon  uncertainty  error  at  opposing  detectors 
are  essentially  un correlated.  The  equivalent  error  in  evaluation  of 

and  R is  then  0.0020°  x ,'2  = 0.0028°  (lo)  or  0.0085°  (3o).  The  output 
error  becomes  0.0424  V for  track  mode  and  0,0106  V for  acquisition  mode. 

The  same  values  can  be  used  for  ?2  and  ^ ' (The  distance  between  look 
points  of  every  other  detector  cell  is  approximately  5383.1  N.Mi.  The 
correlation  of  horizon  uncertainty  at  points  separated  by  this  distance  is 
only  0.12.  It  is  neglected  for  simplicity  of  model.)  For  vehicle  points 
at  nadir  and  spins  at  a yaw  rate  of  W RPM,  the  sensor  detector  scans  the 
earth  surface  at  a rate  of  358.87  W N. Mi/sec.  The  2500  N.  Mi  correlation 
distance  would  then  be  equivalent  to  a correlation  time  of  6.97/W  sec. 

The  sensor  output  error  e(t)  due  to  the  horizon  uncertainty  will  be  modeled 
as  an  exponentially  correlated  noise  as  follows: 

E[e(t)  e(t+T)]  = aQ2  e_T/Tc 


where 


Oq  = 0.0424  V for  track  mode 

= 0.0106  V for  acquisition  mode 

T = 6.97/W  sec 
c 

= 3.485  sec  for  W = 2 RPM 


The  above  model  will  be  used  to  account  for  the  effect  of  horizon  uncertainty 
even  for  cases  where  vehicle  is  not  quite  pointing  at  the  nadir.  (The  change 
of  correlation  time  due  to  the  pointing  offset  from  nadir  by  + 4°  is  expected 
to  be  small.  For  large  offset  from  nadir  the  effect  of  detector  not  seeing 
the  earth  disk  will  be  much  more  significant  than  the  effect  of  horizon 
uncertainty. ) 


In  summary,  the  sensor  error  will  be  simulated  as  an  additive  output  error 
consists  of  a white  noise  component , a bias  component  and  an  exponentially 


correlated  random  noise. 


The  white  noise  is  used  to  simulate  the  effect 


of  detector  and  preanplifier  noises.  The  equivalent  attitude  error  of  the 
white  noise  is  0.0060°  (3a).  The  exponentially  correlated  noise  is  used  to 
simulate  the  effect  of  horizon  uncertainty  error.  The  equivalent  attitude 
error  is  0.0085°  (3a).  The  correlation  time  of  this  error  for  the  vehicle 
spin  rate  of  2 RPM  at  initiation  of  earth  acquisition  is  3.485  sec. 

The  bias  error  is  introduced  to  account  for  the  total  effects  of  electronics 
error,  temperature  gradient  error  and  sensor  alignment  error.  The  equivalent 
attitude  error  is  0.0285°  (3o)  evaluated  as  the  RSS  value  of  each  individual 
e ffects . 

CONCERNS 

o The  detector  geometry  and  output  computations  may  not  be  correct 
since  Barnes  has  two  different  designs.  The  one  used  here  (Figure 
J-5)  is  one  of  their  designs.  The  other  has  detector  pattern  rotated 
22%°  about  the  line  of  sight  from  the  one  shown  here. 

o The  output  update  frequency  for  GPS  is  not  certain.  The  one 
used  here  is  the  one  proposed  by  Barnes  for  DSCS  III. 

0 Earth  image  distortion  and  focus  are  not  modeled.  Excluding  radiance, 
the  earth  image  is  assumed  to  be  a perfect  circle  on  the  focal  surface. 

o The  detector  output  approximation  using  1°  x 1°  squares  causes 
very  little  error  near  null  (normal  operation)  and  up  to  0.15° 
maximum  detector  error  at  extreme  attitudes  (during  acquisition) . 

The  extent  of  the  errors  near  null  and  the  sensor  output  errors 
during  acquisition  has  not  been  analyzed.  The  approximation 
errors  off  null  are  on  the  order  the  errors  due  to  earth  radiance 
variation  and  earth  image  quality.  Smaller  squares  may  be  used  to 
improve  accuracy  but  computer  cost  will  go  up. 
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Figure  J-5.  Detector  Geometry — 2-Axis  Earth  Sensor 
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NUMBERING  OF  X-Y  COORDINATES  FOR  POSITION  OF  1°  x 1°  SQUARES. 
X-Y  COORDINATES  OF  ALL  160  1°  x 1°  SQUARES  ARE  OBTAINED  FROM 
THESE  20  STORED  VALUES 


Figure  *J-5.  Detector  Geometry — 2-Axis  Earth  Sensor  (continued) 
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o Sun  and  moon  interference  are  not  included  in  this  model 


description. 

FLOW  DIAGRAM 

The  combined  earth  sensor  model  flow  diagram  is  shown  in  Figure  J-6.  The 
left  arrow  (*-)  within  boxes  indicates  replacement,  i.e.,  the  expression  on 
the  right  replaces  the  variable  on  the  left.  Mathematical  subscripts  are 
used.  Symbols  used  in  the  flow  diagram  are  defined  in  the  sections  below 
Comments  appear  outside  the  boxes  to  describe  the  operations  within  the 
boxes . 

INITIALIZATION  DATA 

Table  J-2  lists  the  initialization  parameters,  nominal  values,  tolerance, 
and  source  of  data.  This  data  is  to  be  input  at  the  start  of  each  run. 

INPUT  DATA 

Table  j-3  lists  the  input  data  the  earth  sensor  model  expects  from  the 
calling  program  each  time  the  subroutine  is  executed  in  the  normal  mode. 

OUTPUT  DATA 

Table  j-4  lists  the  output  data  from  the  earth  sensor  model  subroutine. 
Other  subroutine  variables  may  be  outputs  if  the  user  desires. 

INTERNAL  VARIABLES 

Internal  variables  used  in  the  flow  diagram  are  defined  in  Table  J-5. 
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Figure  ,T-6.  Earth  Sensor  Flow  Chart  (continued) 
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Figure  J-6.  Earth  Sensor  Flow  Chart  (continued) 
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FigureJ-6.  Earth  Sensor  Flow  Chart  (continued) 
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Figure  J-6.  Earth  Sensor  Flow  Chart  (continued) 
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Compute  area  of  earth 
coverage  on  A and  S 
detectors 


Compute  area  of  earth 
coverage  on  B detectors 


Compute  detactor  output 
voltage  and  save  detector 
coverage . 


Figure  J-6.  Earth  Sensor  Flow 
Chart  (continued) 
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STATE  STATE  + 1 


o* 


Advance  to  the  next  state. 


TYMUP  TYMUP  + 0.125 


0-r 


Increment  time  of  next  update 


TYM  TYM  + At 

I 


Increment  time 


^ RETURN  ^ 


End  of  2-axis  iteration 


(FUNCTION  ^ 
AREA  (XA,  YA)  ) 

_=zz 


FUNCTION  \ b , C , RE  are  common  to  calling 
AREA  (XA,  YA)  ) program.  RD  = 0.0087266  equiv.  radian. 


AREA  ( RE+RD-SQRT  (B-XA)2  + (YA-C)2  ) * — 

7T 


INSERT  (a) 
THRU  (e) 
HERE  FOR 
RADIANCE 
VARIATIONS 


Compute  approx,  area 
of  earth  covering  1°  x 1° 
square . 


Li mi t AREA  to  be 
< 1 and  > 0 


Figure  J-6.  Earth  Sensor  Flow  Chart  (concluded) 


Where 


SIGj  = 0.0020  (deg.,  lo)  1=1,  ....  8 

BSj  - Bias  error  with  worst  case  value  0.0285  (deg.,  3c)  I = 1,  ...,  8. 

REXPj  = Exponentially  correlated  noise  with  standard  deviation 
0.00283°  and  correlation  time  3.485  sec. 

The  function  REXPj  to  be  used  to  generate  the  exponentially  correlated  noise 
is  given  in  the  following: 

FUNCTION  REXPj 

REXPI0  = REXPI 

BETA  = l./TMCNST 

DREXP  = -BETA*REXPI0 

+SQRT  (2. *BETA)*RSIGI*RANDI 

REXPI  = REXPI0  + DREXP 

RETURN 

END 

Where 

TMCNST  = Correlation  Time 
RSIGI  * Standard  Deviation 

RANDI  = Random  nunher  generator  with  unit  variance 
REXPI  to  be  assigned  an  initial  value. 
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TABLE  J-2,  INITIALIZATION  DATA  FOR  EARTH  SENSOR  MODEL 


PARAMETER 

SUBSCRIPT 

UNIT 

NOMINAL 

SOURCE 

TOLERANCE 

SOURCE 

Time  constant  of 

two  axis  earth 

Sec  on  ds 

1.4 

Note  1 

N/A 

None 

sensor  detectors,  T 

Line  of  sight  of 

spin  earth  sensor 

field 

LOS6,  LOS 

1 

-0.9848 

Note  2 

SM 

Not 

speci fied 

2 

-0.9848 

iosY 

1 

0 

2 

0 

'-nV 

1 

0. 1736 

2 

-0.1736 

Transformation  from 

body  to  sensor 

frame:  (TSB) 

Ritw  1 Col.  1 

1 

— 

1.0 

Assumed 

N/A 

None 

1 2 

0 

(nominal  ] 

1 3 

3 

0 

2 l 

0 

2 2 

5 

1.0 

2 3 

6 

0 

3 1 

7 

0 

3 2 

8 

0 

3 3 

9 

1.0 

Calling  t i me 

increment,  t 

Seconds 

L0. 00326 

<0.125 

(spin  model 

(2  axis) 

See 

text 

See 

text 

Center  position 

of  Incremental 

detector  area: 

xi 

1 

Equiv. 

>.25  JO V 

7415  Note 

N/A 

None 

2 

r.idi  an 

(h  2b  3072 

7415  2 A 

1 

0.2705jwn4i  fir.  1 

TABLE  J-2 . INITIALIZATION  DATA  FOR  EARTH  SENSOR  MODEL  (continued) 


PARAMETER  SUBSCRIPT 

UNIT 

NOMINAL 

VALUE 

1 1 

SOURCE  j TOLERANCE  j SOURCE 

4 

0.2705260341 

5 

0.2879  79  3266 

6 

0.2879  793266 

7 

0.305432619 1 

8 

0.3054326191 

9 

0. 3228859116 

in 

0.3228859116 

n 

0.1727787808 

12 

0.160  4 3 7 4 39  3 

13 

0.1851201222 

H 

0.1727787807 

15 

0. 1974614637 

16 

0.  1851201222 

17 

0.2098028052 

18 

0.  1974614637 

19 

0.2221441467 

» 

0.2098028052 

Yi  1 

Equiv. 

0.0087266463 

N/A  j None 

radi  on 

0.0261 799387 

3 

0.0087266463 

0.0261  799387 

5 

0.0087266463 

i 

6 

0.0261799387 

7 

0.0087266463 

8 

0.0261799387 

9 

0.0087266463 

10 

0.0261799387 

11 

0.  1851201222 

12 

0.  1974614637 

13 

0.1974614637 

14 

0.2098028051 

15 

0.2098028051 

16 

, 0.2221441466 

17 

0.2221441466 

18 

0.2344854881 

19 

0.2344854881 

20 





0.2468268296 

1 
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TABLE  J-2.  INITIALIZATION  DATA  FOR  EARTH  SENSOR  MODEL  (concluded) 


PARAMETER 

SUBSCRIPT 

UNIT 

NOMINAL 

VALUE 

SOURCE 

TOLERANCE 

SOURCE 

Noise  Standard 

• • 

Deviation 

SIC  on: 

R1T 

1 

Volts 

0.167 

Note  4 

N/A 

R1A 

2 

Volts 

0.0417 

PIT 

3 

Volts 

0.167 

P1A 

4 

Volts 

0.0417 

R2T 

5 

Volts 

0.167 

R2A 

6 

Volts 

0.0417 

P2T 

7 

Vo  1 ts 

0.167 

P2A 

8 

Volts 

0.0417 

TPIPj 

9 

Sec . 

0 

Not 

Sped  f- 
ied 

Not 

Sped  f- 
led 

tpip2 

10 

Sec. 

0 

Not 

Spec  1 f- 
ied 

Not 

Sped  fr- 
ied 

Two  axis  output 

scale  factor 

Track,  SF 

1 

Volt/sq . 
deg. 

1.25 

Note  3 

0.0625 

Note  3 

Acq . , SF 





Vol  t /sq . 
deg. 

0. 3125 

0.025 

NOTES  TO  TABLE  M-2 
l.  Barne 


Barnes  presentation  to  Honeywell  in  St.  Petersburg  on  3-1  >-75. 
Rockwell  Spec  MC  4 32  -02  1 4,  Rev.  H.,  12-11-74. 

Computed  from  data  in  Rockwell  Spec  MC432-0214,  Rev.  B,  12-11-7 
Track  seal*  factor 


1 degree  attitude  - 4 sq . deg.  coveraRe  * 5 t 0.25V. 

’ * 1.25  t 0.0625  volt/sq.  deg. 

. sq.  deg. 

Acquisition  scale  factor 

4 degree  attitude  - 16  sq.  deg.  coverage 

-M-  - - - 0.3125  - 0.025  «oll/s(.  *■«. 

4 sq  . de  g . 

Confuted  from  data  in  Rockwell  Spec  MC  432-0214,  Rev.  B,  12-11-74 


♦0.1  , 3 s i gma 

0.1  x 5 volt /deg  x 1/3  ■ 0.167  volt,  1 sigma 

Acquisition  mode 

♦0.1°,  3 sigma 

0.1  x 1.25  volt/deg  x 1/3  • 0.0417  volt,  1 sigma 
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TABLE  J-3 


INPUT  DATA  FOR  EARTH  SENSOR  MODEL 


DEFINITION 

FLOW 

DIAGRAM  SYMBOL 

DIMENSION 

TYPE 

UNITS 

Mode  control  flag 

MODE 

In  te  ge  r 

0 = initialize 

1 = spin  opera- 

tion 

2=2  axis 

operation 

Satellite  position 

vector 

R ,R  ,R  ; R 
x y z 

3 

Real 

N Mi 

T rans  formation 

from  body  to 

inertial  frame 

TIB 

9 

Real 

186 


I 


DEFINITION 

FLOW 

DIAGRAM  SYMBOL 

DIMENSION 

TYPE 

UNITS 

Spin  earth  sensor 

output 

TPIP 

2 

Real 

Seconds 

Pitch  attitude 

in  track  mode 

Pitch  attitude 

PT 

Real 

Volts 

in  acquisition 

mode 

PA 

Real 

Volts 

Roll  attitude  in 

track  mode 

RT 

Real 

Volts 

Roll  attitude  in 

acquisition  mode 

RA 

Real 

Volts 

TAELE  J-5.  INTERNAL  VARIABLES  FOR  EARTH  SENSOR  MODEL 


■ 


DEFINITION 

ROW 

DIAGRAM  SYMBOL 

DIMENSION 

TYPE 

UNITS 

Two  axis  rise /decay 

CON  2 

Real 

Volts  per  sq.  deg. 

increment  for  0.25  sec. 

Two  axis  rise /decay 

CON  3 

Real 

Volts  per  sq . deg. 

increment  for  0.125  sec 

Area  of  incremental 

A 

80 

Real 

Square  degree 

detector  covered  by 

. j 

earth 

Time  since  start  of  run 

TYM 

Doub 

Seconds 

prec 

Time  of  next  2-axis 

TYMUP 

Doub 

Sc .onds 

update 

prec 

Satellite  distance  to 

RM 

Real 

N.  M* 

earth  center 

Satellite  position 

RS  , RS  , RS  ; 

X V 2 

3 

Real 

N.  Mi 

vector  in  sensor 

R 

frame 

Radius  of  earth  image 

RL 

Real 

F.quiv.  radian 

on  focal  surface 

Position  of  earth 

B 

Real 

lquiv.  radian 

Image  on  focal 

surfc^e,  X c opponent 

Position  of  earth 

c 

Real 

Equiv.  radian 

image  on  focal 

surface,  Y coirponent 

State  of  sequential 

STATE 

Intege  r 

processing  wi thin 

sensor 

Area  of  detector 

D 

12 

Real 

Square  degree 

Covered  by  earth 

Area  of  detector 

JSAV 

24 

Real 

Square  degree 

coverage  from 

previous  cycle 

Dynamic  Detector 

v 

24 

Real 

Square  Degree 

output 
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TABLE  J-5.  INTERNAL  VARIABLES  FOR  EARTH  SENSOR  MODEL  (concluded) 


DEFINITION 

FLOW 

DIAGRAM  SYMBOL 

DIMENSION 

TYPE 

[ 



UNITS 

Pitch  1 track 

PIT 

Real 

Volts 

atti tude 

Pitch  2 track 

P2T 

Real 

Volts 

atti  tude 

Roll  1 track 

Rjr 

Real 

Volts 

a 1 1 i t ude 

Roll  2 track 

R2T 

Real 

Volts 

at  ti t ude 

Pitch  l ucquisi- 

P1A 

Real 

Volts 

t i or.  at  t L t ude 

Pitch  2 acquis i- 

P2A 

Real 

Vo  I ts 

tion  attitude 

Roll  1 acquisi- 

R1A 

Real 

Volts 

tion  attitude 

Rol  1 2 acq  dsi- 

R2A 

Real  V 

Volts 

tion  attitude 

Index 

i 

I n te  ge  r 

Index 

j 

I n te  ge  r 

Line  of  sight  of  spin 

sensor  field  in 

inertial  frame 

LOS1 

3 

Real 

Normalized  cross 

SANG 

Real 

product  magnitude  of 

K an  d LOS 

Normalized  dot  product 

CANC 

Real 

of  R and  LOS 

Angle  between  -R  and 

ANC1 

Real 

Radian 

tangent  to  earth 

Angle  between  -R  and 

AN  02 

Real 

Radian 

LOS 

Angle  between  LOS 

AN  03 

Real 

Radian 

and  earth  tangent 

Saved  angle  between 

AN  OS  A V 

2 

Real 

Radian 

earth  tangent  and  LOS 

Half  dimension  of 

1 x 1 square 

RD 

Real 

J 

Equiv  radian. 

(a  constant  preset 
to  0.0087266) 

Area  of  earth  cover- 

AREA 

Real 

Square  Degree 

ing  1 x 1 square 

Function 
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APPENDIX  K 

WHEEL-MOMENTUM  COUPLING 


APPENDIX  K 


WHEEL-MOMENTUM  COUPLING 


To  examine  the  effects  of  wheel-momentum  coupling,  we  begin  with  the  rotational 
equations  of  motion,  i.  e.  , 


T = I*  + + V + HW  = TJ  + Td 


where 


f = inertia  tensor 


IITI,  = wheel  momentum 

w 

T = jet  control  torque 
J 


8cu 


(K-l) 


(K-2) 


(K-3) 
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* 


0 


0 


Hw,/I 


+ 0 ID  + 0 

e 


.-"wA 


i/i  o n , + t + t> 

y J W d 

0 1/1 


where 


a;  = earth  rate 
e 


We  have  two  effects  to  examine--one  static  and  the  other  dynamic.  The  static  effect  is 
that  wheel  momentum  induces  wheel  momentum  coupling  torques  in  the  roll  and  yaw 
axes.  If  necessary,  this  error  source  can  be  compensated  for  by  using  tachometer 
feedback  to  generate  wheel  commands  in  the  cross  axis,  i.e.. 


Tw  - -I  0 cu 

Wc  W e 


where 


1„,  = wheel  inertia 
w 


= w^ee^-sPeef*  estimate  from  tachometer 


Assuming  that  wheel  speed  is  known  to  within  one  percent  of  its  true  value,  the  effective 
disturbance  due  to  wheel-momentum  coupling  is  reduced,  i.  e. , 


llWMC 


0 <»  - 

e W 


I,.,  0 cd 

W e 


"^max  ^ WWmax* 


L 


where 


“ nw  = true  minus  estimated  wheel  speed 

The  other  effect  of  wheel-momentum  coupling  is  a dynamic  one.  It  introduces  nutation 
coupling  between  all  axes.  The  maximum  frequency  of  this  nutation  coupling  is  given  by 

0 = [Hw  + (I  - I )(u  ]/l 

xz  Wmax  x y eJ  z 

^xy  ^max^Vyy  (K-7) 

r‘yz  = HWmax/VlzIy 

The  effect  of  this  nutation  coupling  on  control  loops  which  are  designed  to  be  independent 
should  cause  no  problem,  provided  that  these  nutation  frequencies  are  well  below  the 
crossover  frequency  of  the  wheel  loop. 


